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This article reviews recent progress of QCD phase structure, including color superconductor at 
high baryon density and strongly interacting quark-gluon plasma (sQGP) at high temperature cre- 
ated through relativistic heavy ion collision. A brief overview is given on the discovery of sQGP at 
RHIC. The possibility of locating the critical end point (CEP) from the property of bulk viscosity 
over entropy density is discussed. For the phase structure at high baryon density, the status of the 
unconventional color superconducting phase with mismatched pairing is reviewed. The chromomag- 
netic instability, Sarma instability and Higgs instability in the gapless color superconducting phase 
are clarified. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is an asymptoti- 
cally free theory [l| and regarded as the fundamental 
theory of quarks and gluons. At very high energies, in- 
teraction forces become weak, thus perturbation calcu- 
lations can be used. The perturbative QCD predictions 
have been extensively confirmed by experiments, while 
QCD in the non-perturbative regime is still a challenge 
to theorists. The fundamental quarks and gluons of QCD 
have not been seen as free particles, but are always con- 
fined within hadrons. It is still difficult to construct the 
hadrons in terms of nearly massless quarks and gluons. 
The observed baryon spectrum indicates that the (ap- 
proximate) chiral symmetry is spontaneously broken in 
the vacuum. As a result, the eight pseudoscalar mesons 
TT, K and t] are light pseudo-Nambu-Goldstone bosons, 
and the constituent quark obtains dynamical mass, which 
contributes to the baryon mass. Besides conventional 
mesons and baryons, QCD itself does not exclude the ex- 
istence of the non-conventional states such as glueballs, 
hybrid mesons and multi-quark states 

Since 1970s, people have been interested in QCD at ex- 
treme conditions. It is expected that the chiral symme- 
try can be restored, and quarks and gluons will become 
deconfined at high temperatures and/or densities 0-0]. 
Fig. [T]is the typical QCD phase diagram, which shows 
the system is in deconfined quark-gluon plasma phase at 
high temperature, and in color superconducting phase at 
high baryon density. 

Results from lattice show that the quark-gluon plasma 
(QGP) does exist. For the system with zero net baryon 
density, the deconfinement and chiral symmetry restora- 
tion phase transitions happen at the same critical tem- 
perature At asymptotically high temperatures, e.g., 
during the first microseconds of the "Big Bang", the 
many-body system composed of quarks and gluons can be 
regarded as an ideal Fermi and Boson gas. It is believed 
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FIG. 1: QCD phase diagram at finite temperature and 
baryon density. 



that the "little Bang" can be produced at RHIC and 
LHC. Recently, it was shown that the new state of mat- 
ter produced at RHIC is far away from the asymptotically 
hot QGP, but in a strongly coupled regime. This state is 
called strongly coupled quark-gluon plasma (sQGP)B. 
For most recent reviews about QGP, e.g., see Ref. pj|-[l3|. 

Studying QCD at finite baryon density is the tra- 
ditional subject of nuclear physics. The behaviour of 
QCD at finite baryon density and low temperature is 
central for astrophysics to understand the structure of 
compact stars, and conditions near the core of collaps- 
ing stars (supernovae, hypcrnovac). Cold nuclear mat- 
ter, such as in the interior of a Pb nucleus, is at T = 
and iiB — rriN ~ 940MeV. Emerging from this point, 
there is a first-order nuclear liquid-gas phase transition, 
which terminates in a critical endpoint at a temperature 
^ lOMcV [l^. If one squeezes matter further and fur- 
ther, nuleons will overlap. Quarks and gluons in one 
nucleon can feel quarks and gluons in other nucleons. 
Eventually, deconfinement phase transition will happen. 
Unfortunately, at the moment, lattice QCD is facing the 
"sign problem" at nonzero net baryon densities. Our 
understanding at finite baryon densities has to rely on 
effective QCD models. Phenomenological models indi- 
cated that, at nonzero baryon density, the QGP phase 
and the hadron gas are separated by a critical line of 
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roughly a constant energy density e^r — IGeV/fm^ 

In the case of asymptotically high baryon density, the 
system is a color supcrconducor. This was proposed by 
Frautschi [l^ and Barrois [l^. Based on the Bardeen, 
Cooper, and Schrieffer (BCS) theory because there 
is a weak attractive interaction in the color antitriplet 
channel, the system is unstable with respect to the for- 
mation of particle-particle Cooper-pair condensate in the 
momentum space. Detailed numerical calculations of 
color superconducting gaps were firstly carried out by 
Bailin and Love [l^ ■ They concluded that the one-gluon 
exchange induces gaps on the order of 1 MeV at several 
times of nuclear matter density. This small gap has little 
effect on cold dense quark matter, thus the investigation 
of cold quark matter lay dormant for several decades. 
It was only revived recently when it was found that the 
color superconducting gap can be of the order of 100 MeV 
[igj . which is two orders larger than early perturbative 
estimates in Ref. [l^. For this reason, the topic of color 
superconductivity stirred a lot of interest in recent years. 
For review articles on the color superconductivity, see for 
example, Refs. [20j . 

In this article, I will focus on recent progress of sQGP 
created at RHIC and color superconducting phase struc- 
ture at intermediate baryon density regime. The outline 
of this article is as follows: I will give a brief overview 
on the discovery of sQGP at RHIC in Sec. |TT1 Then 
introduce the status of the color superconducting phase 
especially the gapless color superconducting phase in Sec. 
mil At last, I will give a brief outlook in Sec. |IVl 
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FIG. 2: Elliptic flow V2 as a function of pr for different values 
of Ts/tq. The figure taken from Ref. [2g| . 



with the flow velocity, e,p the energy density and 
pressure density, respectively. 

In the Navier-Stokes hydrodynamics, the energy mo- 
mentum tensor decomposes into ideal and dissipative 
parts as 




ideal 



(3) 



with 



II. STRONGLY INTERACTING 
QUARK-GLUON PLASMA(SQGP) 

A. Discovery of sQGP at RHIC 

Studying Quantum chromodynamics (QCD) phase 
transition and properties of hot quark matter at high 
temperature has been the main target of heavy ion col- 
lision experiments at the Relativistic Heavy Ion collider 
(RHIC), the forthcoming Large Hadron Collider (LHC) 
and FAIR at GSI. 

The deconfined quark-gluon plasma, if it can be cre- 
ated through heavy-ion collisions, is an intermediate 
state and cannot be measured directly. In experiment, 
the detector can only measure the freeze-out hadrons. In 
order to extract the property of the intermediate state, 
hydrodynamics is often used to simulate the evolution of 
the fluid. 

The hydrodynamical equations of motion are the local 
conservation laws of energy-momentum and net charge 



'qiV^u" + Vu^ - ^A^'^Vau" + CA^'^V^u". (4) 



(1) 



In ideal hydrodynamics, the energy-momentum tensor 
takes the form of 



ideal 



(e + p)u^vy — pg 



flu 



(2) 



Where A^'^ = .g^'' - uf^u" y = A'^'d^, tjX are the 
shear viscosity and bulk viscosity, respectively. 

It was expected that deconfined quark matter formed 
at high temperature should behave like a gas of weakly 
interacting quark-gluon plasma (wQGP). The perturba- 
tive QCD calculation gives a large shear viscosity in the 
wQGP with r//.s ~ 0.8 for = 0.3 [U. Therefore, 
it turned out as a surprise that the RHIC data of el- 
liptic flow V2 can be described very well by requiring a 
very small shear viscosity over entropy density ratio r]/s 
[13, H^. Lattice QCD calculation confirmed that rj/s for 
the purely gluonic plasma is rather small and in the range 
of 0.1 - 0.2 [13. 

It is now believed that the system created at RHIC is 
a strongly coupled quark-gluon plasma (sQGP) and be- 
haves like a nearly "perfect" fluid The AdS/CFT 
duality gives a lower bound ij/s = l/47r (2^. Therefore, 
it is conjectured that the sQGP created at RHIC might 
be the most perfect fluid observed in nature. 

However, a perfect fluid should have both vanishing 
shear and bulk viscosities. 

The perturbative QCD calculation gives C/* = 0-02al 
for 0.06 < Us < 0.3 [23 • In the hydrodynamic simula- 
tion used to describe the evolution of the fireball created 
at RHIC, the bulk viscosity ( has often been neglected. 



3 





• 




• 


0.8 




0.8 




0.6 




0.6 
0.4 




0.4 


• 


0.2 
0.0 


• 

• 

• . 


0.2 


• 


1.02 1.04 1.06 1.08 1.10 1.12 1 


0.0 


• 

s 






1 




2 3 4 5 



FIG. 3: The bulk viscosity over entropy density ratio as a 
function of scaled temperature T /Tc- The figure is taken from 
Ref.m. 



The zero bulk viscosity is for a conformal equation of 
state and also a reasonable approximation for the weakly 
interacting gas of quarks and gluons. However, recent 
lattice QCD results show that the bulk viscosity over en- 
tropy density ratio C/s rises dramatically up to the order 
of 1.0 near the critical temperature Tc [29l - l3ll |. (There 
are still some subtle issues to determine the bulk viscos- 
ity of QCD through calculating the correlations of the 
energy-momentum tensor on the lattice, see more de- 
tailed discussion in Ref. [32|.) The sharp peak of bulk 
viscosity at has also been observed in the linear sigma 
model [3^ and in the real scalar model [131 ■ The increas- 
ing tendency of C,/ s has been shown in a massless pion 
gas dl] and in the NJL model below [11]. The large 
bulk viscosity near phase transition is related to the non- 
conformal equation of state [13, HI] , and the correlation 
between the bulk viscosity and the conformal anomaly 
has been investigated in Ref. [1^ . 

The sharp rise of the bulk viscosity will lead to the 
breakdown of the hydrodynamic approximation around 
the critical temperature. The effect of large bulk vis- 
cosity on hadronization and freeze-out processes of QGP 
created at heavy ion collisions has been discussed in Refs. 
[i^-dl]. The authors of Ref. [13] pointed out the possi- 
bility that a sharp rise of bulk viscosity near phase tran- 
sition induces an instability in the hydrodynamic flow of 
the plasma, and this mode will blow up and tear the sys- 
tem into droplets. Another scenario is pointed out in Ref. 
[29I l42j that the large bulk viscosity near phase transition 
might induce "soft statistical hadronization" , i.e. the ex- 
pansion of QCD matter close to the phase transition is 
accompanied by the production of many soft partons, 
which may be manifested through both a decrease of the 
average transverse momentum of the resulting particles 
and an increase in the total particle multiplicity. 



B. Searching for the critical end point 

At small baryon chemical potential /x, for QCD with 
two massless quarks, the spontaneously broken chiral 
symmetry is restored at finite temperature, and it is 
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FIG. 4: Searching for the critical end point at RHIC. 



shown from lattice QCD [44| and effective QCD mod- 
els [4^ that this phase transition is of second order and 
belongs to the universality class of 0(4) spin model in 
three dimensions [i^. For real QCD with two quarks of 
small mass, the second order phase transition becomes a 
smooth crossover at finite temperature. At finite baryon 
chemical potential, there are still no reliable results from 
lattice QCD due to the severe fermion sign problem. 
However QCD effective models [45| suggest that the chi- 
ral phase transition at finite ji is of first order. It is ex- 
pected that there exists a critical end point (CEP) in the 
T — n QCD phase diagram. The CEP is defined as the 
end point of the first order phase transition, and belongs 
to the Z{2) Ising universality class |47|. The signature 
of CEP has been suggested in Refs. [4^. The precise 
location of the CEP is still unknown. In the future plan, 
RHIC is going to lower the energy and trying to locate 
the CEP as shown in Fig. |31 

Recently, the authors of Ref. [i^, [13] suggested us- 
ing the shear viscosity over entropy density ratio rj/s to 
locate the CEP by observing the ratio oi rj/s behaves 
differently in systems of water, helium and nitrogen in 
first-, second-order phase transitions, see the system of 
water for example in Fig. [5] The ratio of 77/s shows a 
cusp at Tc for second order phase transition, and a shal- 
low valley near Tc for cross-over, and shows a jump at Tc 
for first-order phase transition. 

Due to the complexity of QCD in the regime of strong 
coupling, results on hot quark matter from lattice cal- 
culation and hydrodynamic simulation are still lack of 
analytic understanding. In recent years, the anti-de Sit- 
tcr/conformal field theory (AdS/CFT) correspondence 
has generated enormous interest in using thermal A/" = 4 
super- Yang-Mills theory (SYM) to understand sQGP. 
The shear viscosity to entropy density ratio 77/s is as 
small as l/47r in the strongly coupled SYM plasma [27j . 
However, a conspicuous shortcoming of this approach is 
the conformality of SYM: the square of the speed of sound 
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FIG. 5: The shear viscosity over entropy density ratio jy/s in 
the water system. The figure is taken from Ref . [50| . 



always equals to 1/3 and the bulk viscosity is always 
zero at all temperatures in this theory. Though C,/satTc 
is non-zero for a class of black hole solutions resembling 
the equation of state of QCD, the magnitude is less than 
0.1 [51|, which is too small comparing with lattice QCD 
results. 

An alternative nonperturbative approach to study 
QCD phase transition is by using effective models. In 
the following, we investigate the thermodynamical and 
transport properties in two toy models, one is the sim- 
plest real scalar model [33,[52|, the other is more relativis- 
tic QCD effective model, i.e, the Polyakov-linear-sigma 
model (PLSM) [H^l, which can describe chiral phase tran- 
sition as well as deconfinement phase transition success- 
fully. 

Real scalar model 

We introduce the real scalar theory including the sextet 
interaction which is described by the Lagrangian 



1 



4 



c(l)^ + H<i). (5) 



When H ^ 0, this theory is invariant under (f> ^cf) and 
has a Z2 symmetry, here a, b, c are model parameters, 
which determine the vacuum properties. The system 
at finite temperature will be evaluated in the Cornwall- 
Jackiw-Tomboulis (CJT) formalism (H^]. We will discuss 
the following four cases: 1) c = 0, 6 > 0, a > 0, iJ = 0, 
the system is always in the symmetric phase. 2) c ~ 
0,b > 0, a < 0, H = 0, the vacuum at T = breaks 
the Z2 symmetry spontaneously, and the symmetry is re- 
stored at higher T with a second-order phase transition. 
3) c = 0, 6 > 0, a < 0, if 7^ 0, the Z{2) symmetry is ex- 
plicitly broken, and the system will experience a crossover 
at high temperature. 4) c > 0, 6 < 0, a > 0, iJ = 0, the 
broken symmetry is restored at high T with a first-order 
phase transition. 

If symmetry is spontaneously broken in the vacuum, 
(j) has a vacuum expectation value 0, (in the case of no 



It is noticed that the new field (j) obtains a tree-level mass 
of m§ = a + 3b (f)'^ + 5c (j)'^. The induced interaction 
terms including the cubic interaction term with coupling 
strength bcp + 10/3c^^, the quartic term with coupling 
strength 6/4-1- 5/2c(p, the quintic term with coupling 
strength ccj), and the six-point interaction term with cou- 
pling strength 1 / 6c. 

Assuming translation invariance, we consider effective 
potential 51 instead of effective action F, these two quan- 
tities are related via: 



F = ^—n, 

T ' 



(8) 



where V is the 3-volume of the system. The effective 
potential in the CJT formalism reads 



\nG-\K) + G-^{K)G{K) - l] (9) 



where no(0) = is the tree- level potential, and 

G'(Go) is the full(tree-level) propagator: 

G-i(K,0) = -A-2+m2(0) ,G~\K,4>) = -K^+mli^) . 

(10) 

In the Hartree approximation, the momentum depen- 
dent contributions are neglected, denotes the contri- 
bution from two-particle irreducible diagrams, and takes 
the form of 



K 



G{K) 
(11) 

The self-consistent one- and two-point Green's func- 
tions satisfy 
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SG 



= 



(12) 



b=d,,G=G 



All thermodynamical information of the system is con- 
tained in the grand canonical potential 57, evaluated at 
the mean field level. The entropy density s is determined 
by taking the derivative of effective potential with respect 
to the temperature, i.e. 



-dn{(j>)/dT. 



(13) 



5 



As the standard treatment in lattice calculation, we in- 
troduce the normalized pressure density p which is nor- 
malized to vanish at T = /i = and the energy density e 

as 



P 



-P 



Ts. 



(14) 



The equation of state p{£) is an important input into 
hydrodynamics. The square of the speed of sound is 
related to p/e and has the form of 



where 



de Tds/dT C^'' 



C,, = de/dT, 



(15) 



(16) 



is the specific heat. At the critical temperature, the 
entropy density as well as the energy density change 
most quickly with temperature, thus one expect that 
should have a minimum at Tc- 

The shear viscosity r] is calculated by using the Boltz- 
mann equation [sHi- The two-particle elastic scattering 
amplitude, which governs particle collisions in the Boltz- 
mann equation, is 
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FIG. 6: The shear viscosity over entropy density 7;/s as a func- 
tion of the temperature T, for cases with a second-order phase 
transition (solid curve), a crossover (dash-dotted curve), and 
with no phase transition for massive field (dashed curve) 
and massless field (dotted curve). The figure is taken from 
Ref.lH. 
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(17) 



where s,t and u are Mandelstam variables, and A3 = 
6Mb+^(l>l + Wc Jj^ G{K, 0)) and A4 = 12(| 5c(/.g 
5c Jj^ G{K, (j))) are effective couplings. 

The shear viscosity over entropy density ratio yy/s in 
the real scalar model in shown in Figl6]and[7]for different 
orders of phase transitions. There is clearly a qualitative 
difference in the rj/s behavior between cases with and 
without a phase transition. It is seen that rj/s shows a 
cusp at Tc for the case of 2nd-order phase transition, a 
shallow valley near Tc for crossover, and shows a jump 
at Tc for the case of Ist-order phase transition. This 
behavior is qualitatively the same as that in the classic 
systems such as the in H2O system as shown in Fig. [5] If 
there is no phase transition, rj/s is always monotonically 
decreasing. 

The bulk viscosity is related to the correlation function 
of the trace of the energy- momentum tensor 



C = — lim — 

9 uj->-a ui 



dt 



(18) 
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FIG. 7: The shear viscosity over entropy density -q/s as a 
function of the temperature T, for the case of 1st order phase 
transition. 



According to the result derived from low energy theorem, 
in the low frequency region, the bulk viscosity takes the 
form of m, 



c 
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9 Wo 
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9 Wo 



T^4^^ + 16|e„| 



dT T4 

{-16e + 9TS + TCy + 16|e„|} . (19) 



with the negative vacuum energy density Sy — Qy ~ 
CI{4')\t=q, and the parameter ujq = tooiT) is a scale at 
which the perturbation theory becomes valid. From the 
above formula, we can see that the bulk viscosity is pro- 
portional to the specific heat Cy near phase transition, 
thus ^/s behaves as l/C^ near Tc in this approximation. 

The bulk viscosity over entropy density ratio C/s as a 
function of T is shown in Figs l8l9lTUl andITT] It is shown 
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FIG. 8; The bulk viscosity over entropy density ("/s as a func- 
tion of T for the case without phase transition in the real 
scalar model. The figure is taken from Ref . [33| . 
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FIG. 11: The bulk viscosity over entropy density ^/s as a 
function of T for a Ist-order phase transition in the real scalar 
model. The figure is taken from Ref.fsJ. 
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FIG. 9: The bulk viscosity over entropy density ("/s as a func- 
tion of T for a 2nd-order phase transition in the real scalar 
model in the real scalar model. The figure is taken from 
Ref.m. 
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FIG. 10: The bulk viscosity over entropy density ("/s as a 
function of T for the case of crossover (the solid line). The 
figure is taken from Ref. [3^ . 



that in the case without symmetry breaking, the bulk 
viscosity over entropy density ^/s decreases monotoni- 
cally with the increase of the temperature. In the case of 
2nd-order phase transition, (^/s decreases with T at low 
temperature region, then rises up at the critical tempera- 
ture Tc and shows an upward cusp, and decreases further 
in the temperature T > T^. In the case of crossover, it is 
observed the cusp behavior of ^/s is washed out. In the 
case of Ist-order phase transition, C/s shows divergent 
behavior at T^. 

From Refs.fsol. [52j . we know that rj/s shows a shal- 
low valley in the case of crossover and a jump at Tc for 
first-order phase transition. But it is hard to distinguish 
whether the system experiences a crossover or first-order 
phase transition just from the value of 77/s extracted from 
the elliptic flow V2- 

From our results in the real scalar model, it is found 
that the ratio of ^/s shows a very sharp peak at Tc in the 
case of first order phase transition, and there is no obvi- 
ous change of C/s for crossover. As pointed out in Ref. 
[40| that a sharp rise of bulk viscosity near phase tran- 
sition induces an instability in the liydrodynamic flow 
of the plasma, and this mode will blow up and tear 
the system into droplets. Therefore, one can distinguish 
whether the system experiences a first order phase tran- 
sition or a crossover from observables at RHIC exper- 
iments. This result supports the idea of using C/s to 
locate the CEP as suggested in Ref. [l^l ■ 

The Polyakov-linear-sigma model 

We have shown the behavior of shear viscosity over 
entropy density t]/s and bulk viscosity over entropy den- 
sity C/s near Tc for different orders of phase transitions 
in a toy model, i.e, the real scalar model. In the follow- 
ing, we use a more realistic QCD effective model, i.e, the 
Polyakov-linear-sigma model (PLSM) , which is described 
by the Lagrangian [56j 



£^Ccl^^ral-U{(j>,(|)*,T) (20) 
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where we have separated the contribution of chiral de- 
grees of freedom and the Polyakov loop. The chiral part 
of the Lagrangian, Cchirai = + Cm consists of the 
fermionic part 

A = Y.^^fin^'D^ - gTaiaa + ij^T^aM'f (21) 

/ 

and the purely mesonic contribution 

Cm = Tr(9^$^9^$-m2$t$)_Ai[Tr($^$)]2 
-A2Tr($t$)2 + c[Dct($) + Det($t)] 
+Tr[i7($ + $t)]^ (22) 

the sum is over the three flavors (f=l,2,3 for u, d, s). 
In the above equation we have introduced a flavor-blind 
Yukawa coupling g of the quarks to the mesons and the 
coupling of the quarks to a background gauge field ~ 
SfioAo via the covariant derivative = — lA^. The 
$ is a complex 3x3 matrix and is defined in terms of 
the scalar cTq and pseudoscalar tTq meson nonets, 



(23) 



The 3x3 matrix H breaks the symmetry explicitly and 
is chosen as 



H — Taha, 



(24) 



where ha are nine external fields. The Ta = \a/'2- are 
the generators of the U{?>) symmetry, Aa are the Gell- 

Mann matrices with Aq = \J^'^- The Ta are normalized 
to Tr(raT,,) = (5q6/2 and obey the [/(3) algebra with 
[Ta,Tiy] = ifabcTc and {Ta.Tb} = datcTc respectively, 
here fabc and dabc for a,b,c = 1,...,8 are the standard an- 
tisymmetric and symmetric structure constants of SU{3) 
group and 



fa 



to 



0, 



dabO — \l '^^ab- 



(25) 



The quantity U{(f>, (p* ,T) is the Polyakov-loop effec- 
tive potential expressed by the dynamics of the traced 
Polyakov loop 



4> - (Tr,i)/iV„ r - (Tr,Lt)/7v^. 



(26) 



The Polyakov loop L is a matrix in color space and ex- 
plicitly given by 



L{x) = Vexp 



13 

I I dTA4{x,T) 




(27) 



with l3 ~ l/T being the inverse of temperature and 
A4 = iA°. In the Polyakov gauge, the Polyakov-loop ma- 
trix can be given as a diagonal representation [s^- The 
coupling between Polyakov loop and quarks is uniquely 
determined by the covariant derivative I?^ in the PLSM 



Lagrangian in Eg. ((20)) . and in the chiral limit, this La- 
grangian is invariant under the chiral flavor group, just 
like the original QCD Lagrangian. The trace of the 
Polyakov-loop, and its conjugate (j)* can be treated as 
classical field variables in this work. 

The temperature dependent effective potential 
U{(j), (j)* ,T) is used to reproduce the thermodynamical 
behavior of the Polyakov loop for the pure gauge 
case in accordance with lattice QCD data, and it has 
the Z{3) center symmetry like the pure gauge QCD 
Lagrangian. In the absence of quarks, we have </> = (/)* 
and the Polyakov loop is taken as an order parameter 
for deconfinement. For low temperatures, U has a 
single minimum at = 0, while at high temperatures 
it develops a second one which turns into the absolute 
minimum above a critical temperature Tq, and the 
Z(3) center symmetry is spontaneously broken. In this 
paper, we will use the potential U (0, </>* , T) proposed in 
Ref.[58|, which has a polynomial expansion in (j> and (p*: 



rp4 



with 



62 (r) = flo + ai 



02 



03 



(28) 



(29) 



A precision fit of the constants , bi is performed to re- 
produce the lattice data for pure gauge theory thermo- 
dynamics and the behavior of the Polyakov loop as a 
function of temperature. The corresponding parameters 
are 

ao ^ 6.75, ai ^ -1.95, 02 = 2.625, 

as = -7.44, 63 = 0.75, 64 = 7.5. (30) 

The critical temperature Tq for deconfinement in the pure 
gauge sector is fixed at 270 MeV, in agreement with the 
lattice results. 

we obtain the thermodynamical potential density as 



with the quarks and antiquarks contribution 

n,-,^^-2TNq [^{ 



(31) 



(27r)3 

ln[l + 3(0 + 0*e-(^«-^)/^)e-(^'-^'/^ + e'^f^^"^)/^] 
+ln[l + 3(0* + 0e-(^'+^)/^)e-(^'+^)/^ + e-3(^J.+M)/T]| 
d'^p 



-2TN, 



ln[l + 3(0 + 0*e-(s.-P)/T)g-(B,-/.)/T ^ g-3(i=;.-^)/T] 
+ln[l + 3(0* + 0e-(^^=+^')/^)e-(^=+^)/^ + e-3(^=+A')/T]|^ 

(32) 
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Here, Nq = 2, = 1, and Eg = ^-j? 

lence quark and antiquark energy for u and d quarks, for 
strange quark s, it is Eg = y7^^ + mf, and m^, is 
the constituent quark mass for it, d and s. The purely 
mesonic potential is 

C/(cr:r,cry) = —{al + al)~hxax~hyay-^^al(jy 

(33) 



?7i^ is the va- 



Minimizing the thermodynamical potential in Eq. (j3ip 
with respective to ay, (j) S'lid 0*, we obtain a set of 
equations of motion 



da.r. 



= 0, 



da„ 



0, 



The set of equations can be solved for the fields as func- 
tions of temperature T and chemical potential /i, and 
the solutions of these coupled equations determine the 
behavior of the chiral order parameter ax, ay and the 
Polyakov loop expectation values 4>, (jf as a function of 
T and /i. 

Fig ll2l shows the pressure density over energy density 
p/s, which is represented in terms of equation-of-state 
(EOS) parameter, at zero density and finite density, re- 
spectively. We observe that the pressure density over en- 
ergy density increases with temperature and saturates at 
high temperature. Both the linear sigma model and the 
Polyakov linear sigma model give very similar results at 
high temperature, the pressure density over energy den- 
sity p/e saturates at a value smaller than 1/3. Another 
common feature of the p/e in the linear sigma model 
and the Polyakov linear sigma model is that there is a 
bump appearing at low temperature region, which is also 
observed in the lattice result. Around the critical tem- 
perature Tc, the pressure density over energy density p/e 
shows a downward cusp. However, the minimum value 
of the p/e around Tc is 0.2 in the linear sigma model, 
which is much larger than the result from the Polyakov 
linear sigma model and the lattice QCD data. For the 
Polyakov linear sigma model, the minimum of p/e around 
Tc is 0.075, which is consistent with the lattice QCD data 

In Fig |131 we plot the bulk viscosity over entropy den- 
sity ratio C/s as a function of the temperature for zero 
chemical potential. It is shown that, at zero chemical po- 
tential ^ = 0, the bulk viscosity over entropy density ^/s 
decreases monotonically with the increase of the temper- 
ature in both the Polyakov linear sigma model and lin- 
ear sigma model, and at high temperature, C/s reaches 
its conformal value 0. In [30j, the bulk viscosity over 
entropy density of the three flavor system is extracted 
from lattice result, which is shown by the square. It is 
observed that C/s in PLSM near phase transition is in 
very good agreement with the lattice result in (30j . i.e, it 
rises sharply near phase transition. 




FIG. 12: The equation-of-state parameter w{T) =p{T)/e{T) 
for ^ = MeV. The Polyakov linear sigma model prediction 
(solid line) and the linear sigma model prediction (dash line) 
are compared with Nf —2+1 lattice QCD data for Nt = 6. 
Lattice data taken from Ref. [S^. The figure is taken from 
Ref.lsl. 



(1=0 MeV 




FIG. 13: The bulk viscosity over entropy density ratio (/s 
as a function of the temperature for = MeV. The solid 
line denotes the Polyakov linear sigma model prediction and 
the dashed line denotes the linear sigma model prediction. 
Lattice data taken from Ref . [SOl] . The figure is taken from 
Ref.lsl. 



In Ref. [3J| , we have investigated the equation of state 
and bulk viscosity in the real scalar model and 0(4) 
model in the case of 2nd order phase transition, crossover 
and 1st order phase transition, and we have found that 
the thermodynamic properties and transport properties 
in these simple models near the critical temperature Tc at 
strong coupling are similar to those of the complex QCD 
system. In a more realistic QCD effective model, i.e, the 
Polyakov linear sigma model |53j , we have systematically 
investigated the thermodynamic properties and bulk vis- 
cosity and found these properties match with lattice data 
very well in the case of zero chemical potential. We fur- 
ther evaluate the chiral phase transitions of u, d and s 
quarks and dcconfincmcnt phase transition at finite tem- 
perature and finite density, and show the T — ii phase 
structure of the Polyakov linear sigma model in Fig|14l 
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n (MeV) T/T^ 



FIG. 14: The T — fi phase diagram in the Polyakov hnear 
sigma model. The figure is taken from Ref . [S^] . 



FIG. 15: The bulk viscosity over entropy density ratio (/s 
in the PLSM for different chemical potentials as functions of 
T/Tc. The figure is taken from Ref. [111. 



For the Polyakov linear sigma model, the result 
shows that the critical end point is around (TetHe) ~ 
(188 MeV, 139.5 MeV), which is close to the lattice re- 
sult {Te,he) = (162 ± 2 MeV,/iE = 120 ± 13 MeV) 
[soj . For the linear sigma model without the Polyakov 
loop, the critical end point is located at (Te,^ie) — 
(92.5MeV,216 MeV). The critical chemical potential fiE 
in PLSM is much lower than that in the PNJL model 
with three quark flavors where the predicted critical end 
point is > 300 MeVHIH. 

The chiral phase transition for the strange quark and 
the deconfinement phase transition in the T— fi plane are 
shown by the dash-dotted line and dotted line, respec- 
tively. It is found that with the increase of chemical po- 
tential, the critical temperature for strange quark to re- 
store chiral symmetry decreases. However, for the decon- 
finement phase transition, with the increase of chemical 
potential, the deconfinement critical temperature keeps 
almost a constant around 220 MeV. It can be seen that 
in the Polyakov linear sigma model, there exists two- 
flavor quarkyonic phase [62| at low density, where the 
u, d quarks restore chiral symmetry but still in confine- 
ment, and three-flavor quarkyonic phase at high density, 
where the u, d, s quarks restore chiral symmetry but still 
in confinement. 

Because the Polyakov-loop in the PLSM is not intro- 
duced dynamically, it is difficult to calculate the trans- 
port properties from the Boltzmann equation. However, 
we can use Ea. p9)) to calculate the bulk viscosity. In 
FiglTSl we plot the bulk viscosity over entropy density 
ratio C/s as a function of the temperature for different 
chemical potentials. It shows (/s as function of the scaled 
temperature T/Tc for different chemical potentials with 
/i = 0,80,139.5,160 MeV. We can see that when the 
chemical potential increases up to /j, = 80 MeV, there is 
an upward cusp appearing in(/s right at the critical tem- 
perature Tc- With the increase of the chemical potential, 
the upward cusp becomes sharper, and the height of the 
cusp increases. At the critical end point fiE and when 
fi > fiE for the first order phase transition, C/s becomes 



divergent at the critical temperature. 

As discussed earlier, the sharp rise of the bulk viscosity 
will lead to the breakdown of the hydrodynamic approx- 
imation around the critical temperature, and will affect 
the hadronization and freeze-out processes of QGP cre- 
ated at heavy ion collisions Refs. (40l - |43| . For example, 
in Ref. , it is pointed out that a sharp rise of bulk vis- 
cosity near phase transition might induce an instability 
in the hydrodynamic flow of the plasma, and this mode 
will blow up and tear the system into droplets. Another 
scenario is pointed out in Ref. [H, |42| that the large bulk 
viscosity near phase transition might induce "soft statis- 
tical hadronization" , i.e. the expansion of QCD mat- 
ter close to the phase transition is accompanied by the 
production of many soft partons, which may be mani- 
fested through both a decrease of the average transverse 
momentum of the resulting particles and an increase in 
the total particle multiplicity. Therefore the critical end 
point might be located through the observables which 
are sensitive to the ratio of bulk viscosity over entropy 
density. 

However, it is noticed that the results of bulk viscosity 
in this paper are based on Eq. (jl9p , where the ansatz for 
the spectral function 

^ = (35) 
cj TT 7r(a;2 + u^) ^ ' 

has been used in the small frequency, and wq is a scale 
at which the perturbation theory becomes valid. In our 
calculation, loq = I GeV, its magnitude at Tc is in agree- 
ment with that obtained in ChPT for massive pion gas 
system in Ref. [l^. Qualitatively, the bulk viscosity 
corresponds to nonconformality, thus it is reasonable to 
observe a sharp rising of bulk viscosity near phase tran- 
sition. Ref. [33] has investigated the correlation between 
the bulk viscosity and conformal breaking, and supports 
the results in Ref. [l^, [sO] . The sharp rising of bulk vis- 
cosity has also been observed by another lattice result 
[3l| and in the linear sigma model [s^. However, till 
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now, no full calculation has been done for the bulk vis- 
cosity. The frequency dependence of the spectral density 
has been analyzed in Refs. [6^ and [s^l and the limita- 
tion of the ansatz Eg. ([55]) has been discussed. From Eq. 
(fT9)) . we see that the bulk viscosity is dominated by Cy at 
Tc- If Cv diverges at Tc, the bulk viscosity should also be 
divergent at the critical point and behave as How- 



ever, the detailed analysis in the Ising model in Ref. [64 1 



shows a very different divergent behavior C, ^ t~ 
with z ~ 3 the dynamic critical exponent and u ~ 0.630 
the critical exponent in the Ising system. More careful 
calculation on the bulk viscosity is needed in the future. 



III. COLOR SUPERCONDUCTING PHASES 



Cooper pairing of the electrons breaks the electromag- 
netic gauge symmetry, and the photon obtains an effec- 
tive mass. This indicates the Meissner effect [6^, i.e., a 
superconductor expels the magnetic fields. 

In QCD case at asymptotically high baryon density, 
the dominant interaction between two quarks is due to 
the onc-gluon exchange. This naturally provides an at- 
tractive interaction between two quarks. The scattering 
amplitude for single-gluon exchange in an SU (Nc) gauge 
theory is proportional to 



{Ta)ki{Ta)lj 



{SjkSil - StkSji) (37) 



Nc-1 



{SjkSxi + StkSji). 



It is known that the deconfined cold dense quark mat- 
ter is in color superconducting phase. 

Let us start with the system of free fermion gas. 
Fermions obey the Pauli exclusion principle, which means 
no two identical fermions can occupy the same quantum 
state. The energy distribution for fermions (with mass 
to) has the form of 



1 



, /3 = 1/T, 



(36) 



fi is the chemical potential and 



here Ep = \r[^ + 
T is the temperature. At zero temperature, f{Ep) = 
9{fi — Ep). The ground state of the free fermion gas is 
a filled Fermi sea, i.e., all states with the momenta less 
than the Fermi momentum pf = ^//i^ — to^ are occu- 
pied, and the states with the momenta greater than the 
Fermi momentum pp are empty. Adding or removing a 
single fermion costs no free energy at the Fermi surface. 

For the degenerate Fermi gas, the only relevant fermion 
degrees of freedom are those near the Fermi surface. Con- 
sidering two fermions near the Fermi surface, if there is a 
net attraction between them, it turns out that they can 
form a bound state, i.e.. Cooper pair (65| . The binding 
energy of the Cooper pair A{K) {K the total momen- 
tum of the pair), is very sensitive to K, being a maxium 
where K = 0. There is an infinite degeneracy among 
pairs of fermions with equal and opposite momenta at 
the Fermi surface. Because Cooper pairs are composite 
bosons, they will occupy the same lowest energy quantum 
state at zero temperature and produce a Bose-Einstein 
condensation. Thus the ground state of the Fermi sys- 
tem with a weak attractive interaction is a complicated 
coherent state of p article-particle Cooper pairs near the 
Fermi surface [13|. Exciting a quasiparticle and a hole 
which interact with the condensate requires at least the 
energy of 2 A. 

In QED case in condensed matter, the interaction be- 
tween two electrons by exchanging a photon is repul- 
sive. The attractive interaction to form electron-electron 
Cooper pairs is by exchanging a phonon, which is a col- 
lective excitation of the positive ion background. The 



Where Ta is the generator of the gauge group, and i, j and 
fc, / are the fundamental colors of the two quarks in the 
incoming and outgoing channels, respectively. Under the 
exchange of the color indices of either the incoming or the 
outgoing quarks, the first term is antisymmetric, while 
the second term is symmetric. For Nc = 3, Eq. psp rep- 
resents that the tensor product of two fundamental colors 
decomposes into an (antisymmetric) color antitriplet and 
a (symmetric) color sextet. 



[31^ = [31: 



[61 ' 



(38) 



In Eq. psp . the minus sign in front of the antisymmet- 
ric contribution indicates that the interaction in this an- 
titriplet channel is attractive, while the interaction in the 
symmetric sextet channel is repulsive. 

For cold dense quark matter, the attractive interac- 
tion in the color antitriplet channel induces the conden- 
sate of the quark-quark Cooper pairs, and the ground 
state is called the "color superconductivity". Since the 
diquark cannot be color singlet, the diquark conden- 
sate breaks the local color SU{3)c symmetry, and the 
gauge bosons connected with the broken generators ob- 
tain masses. Comparing with the Higgs mechanism of 
dynamical gauge symmetry breaking in the Standard 
Model, here the diquark Cooper pair can be regarded as 
a composite Higgs particle. The calculation of the energy 
gap and the critical temperature from the first principles 
has been derived systematically in Refs. j67H74| . 

In reality, we are more interested in cold dense quark 
matter at moderate baryon density regime, i.e., ~ 
500MeV, which may exist in the interior of neutron stars. 
It is likely that cold dense quark droplet might be created 
in the laboratory through heavy ion collisions in GSI- 
SPS energy scale. At these densities, an extrapolation of 
the asymptotic arguments becomes unreliable, we have to 
rely on effective models. Calculations in the framework of 
pointlike four-fermion interactions based on the instanton 
vertex as well as in the Nambu-Jona-Lasinio 

(NJL) model 78l-l82| show that color superconductivity 



does occur at moderate densities, and the magnitude of 
diquark gap is around 100 MeV. 
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A. Different color superconducting phases 

Even though the antisymmetry in the attractive chan- 
nel signifies that only quarks with different colors can 
form Cooper pairs, color superconductivity has very rich 
phase structure because of its flavor, spin and other de- 
grees of freedom. In the following, I list some of the 
known color superconducting phases. 

The 2SC phase 

Firstly we consider a system with only massless u and d 
quarks, assuming that the strange quark is much heavier 
than the up and down quarks. The color superconducting 
phase with only two flavors is normally called the 2SC 
phase. 

Renormalization group arguments [13, [13, H^l suggest 
that possible quark pairs always condense in the s— wave. 
This means that the spin wave function of the pair is 
anti-symmetric. Since the diquark condenses in the color 
antitriplet 3c channel, the color wave function of the pair 
is also anti-symmetric. The Pauli principle requires that 
the total wave function of the Cooper pair has to be anti- 
symmetric under the exchange of the two quarks forming 
the pair. Thus the flavor wave function has to be anti- 
symmetric, too. This determines the structure of the 
order parameter 

Af/ = Ae,Je"'^^ (39) 

where color indices a, (3 € {r,g,b) and flavor indices 
i,j € {u,d). From the order parameter Eq. pop. we 
can see that the condensate picks a color direction (here 
the blue direction, which is arbitrarily selected). The 
ground state is invariant under an SU{2)c subgroup of 
the color rotations that mixes the red and green colors, 
but the blue quarks are singled out as different. Thus 
the color SU{3)c is broken down to its subgroup SU{2)c, 
and five of the gluons obtain masses, which indicates the 
Meissner effect j8a |. 

In the 2SC phase, the Cooper pairs are ud — du sin- 
glets and the global flavor symmetry SU{2)l ® SU{2)fi 
is intact, i.e., the chiral symmetry is not broken. There is 
also an unbroken global symmetry which plays the role 
of U{1)b- Thus no global symmetry arc broken in the 
2SC phase. 

The CFL phase 

In the case when the chemical potential is much larger 
than the strange quark mass, we can assume ttIu — ma = 
vris = 0, and there are three degenerate massless flavors 
in the system. The spin-0 order parameter should be 
color and flavor anti-symmetric, which has the form of 

^S' = ^E^'^^^"''' (40) 

where color indices a, /3 e (r, g, h) and flavor indices i, j € 
(m, d, s). Writing e.y/e"'^^ = 8f5^j — SjS^ , we can see 



that the order parameter 

Af = A(<5f <5f - 6f6^) (41) 

describes the color-flavor locked (CFL) phase p ro p osed in 
Ref . nil . Many other different treatments [87H89| agreed 
that a condensate of the form PO]) is the dominant con- 
densate in three-flavor QCD. 

In the CFL phase, all quark colors and flavors partici- 
pate in the pairing. The color gauge group is compl etely 
broken, and all eight gluons become massive [80, |90| . 
which ensures that there are no infrared divergences as- 
sociated with gluon propagators. Electromagnetism is no 
longer a separate symmetry, but corresponds to gauging 
one of the flavor generators. A rotated electromagnetism 
("Q") remains unbroken. 

Two global symmetries, the chiral symmetry and the 
baryon number, are broken in the CFL phase, too. 
In zero-density QCD, the spontaneous breaking of chi- 
ral symmetry is due to the condensation of left-handed 
quarks with right-handed quarks. Here, at high baryon 
density, the chiral symmetry breaking occurs due to a 
rather different mechanism: locking of the flavor rota- 
tions to color. In the CFL phase, there is only pairing 
of left-handed quarks with left-handed quarks, and right- 
handed quarks with right-handed quarks, i.e., 

< - < i^mK, > ■ (42) 

Where L, R indicate left- and right-handed, respectively, 
a,/3 are color indices and i,j are flavor indices. A gauge 
invariant form [9ll . [9^ 

^ A2e„-„e«" (43) 

captures the chiral symmetry breaking. The spectrum of 
excitations in the CFL phase contains an octet of Gold- 
stone bosons associated with the chiral symmetry break- 
ing. This looks remarkably like those at low density. In 
the excitation spectrum of the CFL phase, there is an- 
other singlet U{1) Goldstone boson related to the baryon 
number symmetry breaking, which can be described us- 
ing the order parameter 

< udsuds >^< AA > . (44) 

In QCD with three degenerate light flavors, the spectrum 
in the CFL phase looks similar to that in the hyper- 
nuclear phase at low-density. It is suggested that the low 
density hypcr-nuclcar phase and the hi gh d ensity quark 
phase might be continuously connected [93|. 

Spin-1 color superconductivity 

In the case of only one-flavor quark system, due to the 
antisymmetry in the color space, the Pauli principle re- 
quires that the Cooper pair has to occur in a symmetric 
spin channel. Therefore, in the simplest case, the Cooper 
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pairs carry total spin one. Spin-1 color superconductiv- 
ity was firstly studied in Ref. [3, for more recent and 
detailed discussions about the spin-1 gap, its critical tem- 
perature and Meissner effect, see Refs. [sl [toI. [oi - IOTj . 
For a review, see Ref. [98| . 

Pairing with mismatch: LOFF, CFL-K, g2SC 
and gCFL 

To form the Cooper pair, the ideal case is when the 
two pairing quarks have the same Fermi momenta, i.e., 

PF,i = Pfj withpi?_i — ^ IJ-'p i — rri^, like in the ideal 2SC, 
CFL, and spin-1 color superconducting phases. However, 
in reality, the nonzero strange quark mass or the require- 
ment of charge neutrality induces a mismatch between 
the Fermi momenta of the two pairing quarks. When the 
mismatch is very small, it has little effect on the Cooper 
pairing. While if the mismatch is very large, the Cooper 
pair will be destroyed. The most interesting situation 
happens when the mismatch is neither very small nor 
very large. 

LOFF: In the regime just on the edge of decoupling 
of the two pairing quarks (due to the nonzero strange 
quark mass for the qs Cooper pair with q € (it, d) or 
the chemical potential difference for the ud Cooper pair), 
a "LOFF" (Larkin-Ovchinnikov-Fulde-Ferrell) state may 
be formed. The LOFF state was firstly investigated in 
the context of electron supercon ductivity in the presence 
of magnetic impurities (99l . llOOf . It was found that near 
the unpairing transition, it is favorable to form a state in 
which the Cooper pairs have nonzero momentum. This 
is favored because it gives rise to a regime of phase space 
where each of the two quarks in a pair can be close to its 
Fermi surface, and such pairs can be created at low cost 
in free energy. This sort of condensates spontaneously 
break translational and rotational invariance, leading to 
gaps which vary periodically in a crystalline pattern. The 
crystalline color superconductivity h as been in vestigated 
in a series of papers, e.g., see Refs. [ioMqI. 

CFL-K: The strange quark mass rris induces an effec- 
tive chemical potential fis = Ws/(2pf), and the effects 
of the strange quark mass can be quite dramatic. In the 
CFL phase, the and K° modes may be unstable for 
large value s of the s trange quark mass to form a kaon con- 
densati on Il09l- 112j . In the framework of effective theory 
(oll lo^lTTS - llq .The masses of the Goldstone bosons can 



be determined as 



m^± = =F- 



2pf 



2pF 



2pf 



4A 

J TT 
J TT 



-jjimu + ms)md 



-jjims + md)mu 



1/2 



1/2 



1/2 



(45) 



masses are substantially affected by the strange quark 
mass, the masses of K~ and are pushed up while 
and are lowered. As a result, the and K° become 
massless if ms\crit = 3.03 rri^J^A'^^^. For larger values of 
771s the kaon modes are unstable, signaling the formation 
of a kaon condensate. Recently, it was found t hat i n the 
CFL phase, there also may exist rj condensate [ll7l |. 

g2SC and gCFL: When the /3-equilibrium and the 
charge neutrality condition are required for the two-flavor 
quark system, the Fermi surfaces of the pairing u quark 
and d quark differ by ^e, here fig is the chemical potential 
for electrons. It was found that when the gap parameter 
A < /ie/2, the system will be i n a n ew ground state called 
the gapless 2SC (g2SC) phase [ill. The g2SC phase has 
very unusual tempe ratur e properties |ll9l | and chromo- 
magnetic properties [l20j | . This phase will be introduced 
in more detail in Sec. IIII Bl 

Similarly, for a charge neutral 3-flavor system with a 
nonzero strange quark mass 777^, with increasing 777^, the 
CFL phase transfe rs to a new gapless CFL (gCFL) phase 
when 777^ //i ~ 2 A [l2ll | . The finite temperature property 
of the charge neu tral three -flavor quark matter was in- 
vestigated in Ref. p^ill2^ . Recently, it was shown that 



the kaon condensate shifts the critical strange quark mass 
to h igher values for the appearance of the gCFL phase 
fTH. 



with A = 3A2/(47r2) [Uliil]. It was found that the kaon 



B. Gapless color superconductor 

In this section, I would like to focus on unconventional 
color superconductor with mismatched pairing by taking 
charge neutral two-flavor system as an example. 

It is very likely that the color superconducting phase 
may exist in the core of compact stars, where bulk mat- 
ter should satisfy the charge neutrality condition. This 
is because bulk matter inside the neutron star is bound 
by the gravitation force, which is much weaker than the 
electromagnetic and the strong color forces. Any elec- 
tric charges or color charges will forbid the formation of 
bulk matter. In addition, matter inside neutron star also 
needs to satisfy the /3-equilibrium. 

In the ideal two-flavor color superconducting (2SC) 
phase, the pairing u and d quarks have the same Fermi 
momenta. Because u quark carries electric charge 2/3, 
and d quark carries electric charge —1/3, it is easy to 
check that quark matter in the ideal 2SC phase is pos- 
itively charged. To satisfy the electric charge neutrality 
condition, roughly speaking, twice as many d quarks as 
u quarks are needed. This induces a large difference be- 
tween the Fermi surfaces of the two pairing quarks, i.e., 
fJ'd ~ fJ-u = ~ m/4, where 11, fie are chemical poten- 
tials for quarks and electrons, respectively. Naively, one 
would expect that the requirement of the charge neutral- 
ity condition will destroy the ud Cooper pairing in the 
2SC phase. 

Indeed, the interest in the charge neutral 2SC phase 
was stirred by the paper Ref. |126| . It was claimed in this 
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paper that there wih be no 2SC phase inside neutron star 
under the requirement of the charge neutrahty condition. 
In fact, the authors meant that for a charge neutral three 
flavor system, the 2SC+S phase is not favorable compared 
to the CFL phase. This is a natural result under the 
assumption of a small strange quark mass, even without 
the requirement of the charge neutrality condition. In 
the framework of the bag model, in which the strange 
quark mass is very small, the CFL phase is always the 
ground state for cold dense quark matter, and there is no 
space for the existence of two-flavor quark matter. 

However, there is another scenario about the hadron- 
quark phase transition in the framework of the SU(3) 
NJL model. In the vacuum, quarks obtain their dy- 
namical masses induced by the chiral condensate, u, d 
quarks have constituent mass around 330MeV, while the 
s quark has heavier constituent mass, which is around 
SOOMeV. With the increasing of the bayron density, the 
constituent quark mass starts to decrease when the chem- 
ical potential becomes larger than its vacuum constituent 
mass. In this scenario, s quark restores chiral symmetry 
at a lar ger critical chemical potential than that of u, d 
quarks. If the deconfinement phase transition happens 
sequentially, there will exist some baryon density regime 
for only u, d quark matter and s quark is still too heavy 
to appear in the system. 

It is worth to mention that the effect of the electric 
charge neutrality condition on a three-flavor quark sys- 
tem is very different from that on a two-flavor quark 
system. Because s quark carries —1/3 electric charges, 
it is much easier to neutralize the electric charges in a 
three-flavor quark system than that in a two-flavor quark 
system. However, the color charge neutrality condition 
is nontrivial in a three-flavor quark system, when the 
strange quark mass is not very small. For a detailed con- 
sideration of the ch arge neutral three-flavor system, see 
recent papers Refs. [I2ll4l24 |. 

In the following, we focus on the charge neutral two 
flavor quark system. Motivated by the sequ ential decon- 
flnement scenario, the authors of Ref. [l27l | investigated 
charge neutral quark matter based on the SU(3) NJL 
mod el. T o large extent, their results agree with those in 
Ref. |l26{ . i.e., the CFL phase is more favorable than the 
2SC-I-S phase in charge neutral three-flavor cold dense 
quark matter, and they did not flnd the charge neutral 
2SC phase. 

However, it was found in Ref. |l28| that a charge neu- 
tral two- flavor color superconductin g (N2SC) phase does 
exist, which was conflrmed in Refs. jl29l ll3Cll |. Compar- 
ing with the ideal 2SC phase, the N2SC phase found in 
Ref. [128] has a largely reduced diquark gap parameter, 
and the pairing quarks have different number den sities . 
The latter contradicts the paring ansatz in Ref. |l3l| . 
Therefore, one coul d sug gest that t his p hase is an un- 
stable Sarma state jl32l |. In Ref. |ll8l |. it was shown 
that the N2SC phase is a stable state under the restric- 
tion of the charge neutrality condition. As a by-product, 
which comes out as a very important feature, it was 



found that the quasi-particle spectrum has zero-energy 
excitation in this charge neutral two-flaovr color super- 
conducting phase. Thus this phase is named the "gapless 
2SC(g2SC)" phase. 

The 2-flavor system can be described by the gauged 
Nambu-Jona-Lasinio (gNJL) model, the Lagrangian den- 
sity has the form of 
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with Df^i = df^i — igA^jT"'. Here are gluon flelds, 
= A°/2 are the generators of SU{2>)c gauge group 
with a = 1, • • • ,8. In the gNJL model, the gauge flelds 
are external fields and do not contribute to the dynamics 
of the system. The property of the color superconduct- 
ing phase characterized by the diquark gap parameter is 
determined by the nonperturbative gluon fields, which 
has been simply replaced by the four-fermion interaction 
in the NJL model. Gs and Gd are the quark- antiquark 
coupling constant and the diquark coupling constant, re- 
spectively, q^ = G(f , (f = q^C are charge-conjugate 
spinors, G = 17^7'' is the charge conjugation matrix (the 
superscript T denotes the transposition operation). The 
quark field q = qia with i = u,d and a = r, 5, & is a flavor 
doublet and color triplet, as well as a four-component 
Dirac spinor, r = (r^,T^,r'^) are Pauli matrices in the 



flavor space, where is antisymmetric, and (e)**^ = c , 
(e^)"^ = e"^^ are totally antisymmetric tensors in the 
flavor and color spaces. 

In /^-equilibrium, the matrix of chemical potentials in 
the color-flavor space p. is given in terms of the quark 
chemical potential fi, the chemical potential for the elec- 
trical charge fie and the color chemical potential /ig, 
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The total thermodynamic potential for w, d quarks in 
-equilibrium with electrons takes the form |ll8l . Ill9l 
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where the electron mass was taken to be zero, which is 
sufficient for the purposes of the current study. The sum 
over A runs over all (6 quark and 6 antiquark) quasi- 
particles. The explicit dispersion relations and the de- 
generacy factors of the quasi-particles read 



= E{p)±fj.db, 
= E^{p)±Sfi. 
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FIG. 16: Dispersio n rel ation of the gapless phase. The figure 
is taken from Ref. UJM- 



Gapless excitation in quasi-particle spectrum 

It is instructive to start with the excitation spectrum 
in the case of the ideal 2SC phase when 5^ — 0. With the 
conventional choice of the gap pointing in the anti-blue 
direction in the color space, the blue quarks arc not af- 
fected by the pairing dynamics, and the other four quarsi- 
particle excitations are linear superpositions of Ur,g and 
dr^g quarks and holes. The quasi-particle is nearly iden- 
tical with a quark at large momenta and with a hole 
at small momenta. We represent the quasi-particle in 
the form of Q{quark, hole)^ then the four quasi-particles 
can be represented explicitly as Q{ur,dg), Q(ug,dr), 
Q{dr,Ug) and Q{dg,Ur). When S^j, = 0, the four quasi- 
particles are degenerate, and have a common gap A. If 
there is a small mismatch (Sfi < A) between the Fermi 
surfaces of the pairing u and d quarks, the excitation 
spectrum will change. It is found that dfi induces two 
different dispersion relations, the quasi-particle Q{dg, Ur) 
has a smaller energy gap A — Sfi, and the quasi-particle 
Q{ur, dg) has a larger energy gap A -I- 5//. This is similar 
to the case when the mismatch is in duced by the mass 
difference of the pairing quarks [l33| . 

If the mismatch 5^ is larger than the gap parame- 
ter A, the lower dispersion relation for the quasi-particle 



Q{dg 



will cross the zero-energy axis, as shown in the 



Ur) vanishes at two values of momenta p = fi and 
+ where fi^ = fl± (Sfi)^ — A^. Thus this phase 



right panel of Fig. 1161 The energy of the quasi-particle 

Qidg 
P = H 

is called the gapless 2SC (g2SC) phase. 

Thermal stable charge neutral g2SC state 

An uns table gapless CFL phase has been found in 
Ref. jl33l |. and a similar stable gapless color supercon- 



ductivity could also appear in a cold atomic gas [134 
or in ti, s or d, s q uark matter when the number densi- 
ties are kept fixed jl35| . Also, some gapless phases may 
appe ar du e to P-wave interactions in the cold atomic sys- 
tem jl36l |. However, the gapless 2SC phase is a thermal 
stable state under the charge neutrality condition. 

If a macroscopic chunk of quark matter exists inside 
compact stars, it must be neutral with respect to electric 
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FIG. 17: The eflective potential as a function of the diquark 
gap A calculated at a fixed value of the electric chemical po- 
tential fie = 148 MeV (dashed fine), and the efi'ective poten- 
tial defined along the neutrality line (solid line). The results 
are plotted for — 400 MeV with rj — 0.75. The figure is 
taken from Ref. IllSll. 



as well as color charges. Now, we discuss the role of 
the electric charge neutrality condition. If a macroscopic 
chunk of quark matter has nonzero net electric charge 
density ng, the total thermodynamic potential for the 
system should be given by 



n = n 



Coulomb 



(52) 



where fi, 



Coulomb 



UgV'^^^ {V is the volume of the 



system) is induced by the repulsive Coulomb interac- 
tion. The energy density grows with increasing the vol- 
ume of the system, as a result, it is almost impossible 
for matter inside stars to remain charged over macro- 
scopic distances. So bulk quark matter must satisfy elec- 
tric neutrality condition with flcouiomb\nQ=a ~ 0, and 
^u,d,e\nQ=o IS On the neutrality line. Under the charge 
neutrality condition, the total thermodynamic potential 
of the system is il.\nQ=o = ^u,d,e\nQ=o- 

Here, we want to emphasize that: The proper way to 
find the ground state of homogeneous neutral u, d quark 
matter is to minimize the thermodynamic potential along 
the neutrality line r2|„Q=o = ^u,d,e\nQ=o- This is differ- 
ent from that in the flavor asymmetric quark system, 
where /^-equilibrium is required but is a free param- 
eter, and the ground state for flavor asymmetric quark 
matter is determined by minimizing the thermodynamic 
potential ftu,d,e- At a fixed /ig = 148 MeV and with color 
charge neutrality, the thermodynamic potential is shown 
as a function of the diquark gap by the dashed line in 
Fig. [T71 The minimum gives the ground state of the fla- 
vor asymmetric system, and the corresponding diquark 
gap is A = 0, but this state has negative electric charge 
density, and cannot exist in the interior of compact stars. 



C. Chromomagnetic instability in the g2SC phase 

As we know, one of the most important properties 
of the ordinary superconductor is its Mcissner effect. 
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i.e., a superconductor expels the magnetic field, which 
was discovered by Meissner and Ochsenfeld in 1933 j66j . 
From theoretical point of view, the Messiner effect can 
be explained using the linear response theory. The in- 
duced current j-"'' is related to the magnetic field Aj by 
jznd _ HyA-' , where the response function Hij is the pho- 
ton polarization tensor. The response function ha s tw o 
components, diamagnctic and paramagnetic part [137 1. 
In the static and long-wavelength limit, for the normal 
metal, the paramagnetic component cancels exactly the 
diamagnctic component. While in the superconducting 
phase, the paramagnetic component is quenched by the 
energy gap and producing a net diamagnctic response. 
Thus the ordinary superconductor is a perfect diamag- 
net. 

In color superconducting phases, the gluon self-energy 
(the response function to an external color field), has 
been investigated in the ideal 2SC phase [s^ and in the 
CFL phase [901 ■ The results show that the gauge bosons 
connected with the broken generators obtain masses in 
these phases, which indicate the Meissner screening effect 
in these phases. 

It is very interesting to know the chromomagnetic 
property in the g2SC phase. We studied the g2SC phase 
in the framework of the SU(2) NJL model, and the NJL 
model lacks gluons. As reflection of this, it possesses 
the global instead of gauged color symmetry. In addi- 
tion, there appear five Nambu-Goldstone (NG) bosons 
in the ground state of the model when the color symme- 
try is broken. In QCD, there is no room for such NG 
bosons. However, the NJL model can be thought of as 
the low energy theory of QCD in which the gluons, as 
independent degrees of freedom, are integrated out. The 
gluons could be reintroduced back by gauging the color 
symmetry in the Lagrangian density of the NJL model, 
providing a semirigorous framework for studying the ef- 
fect of the Cooper pairing on the physical properties of 
gluons. 

The existence of the g2SC phase can be regarded as a 
physical and model independent result under the restric- 
tion of local charge neutrality conditio n, th e order param- 
eter for this phase is A < S/j.. In Ref. |l2d |. we calculated 
the gluon self-energy in the g2SC phase. It is found that, 
in this phase, the symmetry broken gauge bosons have 
imaginary Meissner screening masses, which is induced 
by the dominant paramagnetic contribution to the gluon 
self-energy. In condensed matter, this pheno men on is 
called the paramagnetic Meissner effect (PME) [l38{ . and 
has been observed in some high temperature supercon- 
ductors and small superconductors. 

Unavoidably, the imaginary Meissner screening mass 
indicates a chromomagnetic instability of the g2SC 
phase. There are many proposals on resolving the chro- 
momagnetic instability jl39l4l43 [: One is through a gluon 



condensate to stabilize the system, which may not change 
the structure of the g2SC phase. It is also possible 
that the instability drives a new stable ground state, 
which may have a rotational symmetry breaking like in 
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FIG. 18: Squared values of the gluon Debye (upper panel) 
and Meissner (lower panel) screening masses, devided by 
nig — Aasfi^ /3tv, as functions of the dimensionless parameter 
A/S^. The red solid line denotes the results for the gluons 
with A = 1, 2, 3, the green long-dashed line denotes the results 
for the gluons with A = 4, 5, 6, 7, and the blue short-dashed 
line denotes the re sults for the gluon with A = 8. The figure 
is taken from Ref. \l2d\. 



Refs. Il45l |. or even havean inhomogeneous phase 
structure, like a crystal [99l - ll0l| or a vortex [l46| struc- 
ture. 



D. Sarma instability and Higgs instability 

In order to understand the chromomagnetic instability 
in the gapless phases, we extend t he ga uged NJL model 
beyond mean- field approximation jl42l |. The supercon- 
ducting state is characterized by the order parameter 
A(a;), which is a complex scalar field and has the form of 
A{x) = |A(a;)|e''^(^), with |A| the amplitude and ip the 
phase of the gap order parameter. For a homogeneous 
condensate, A(a;) is a spatial constant. The fluctuations 
of the phase give rise to the pseudo Nambu-Goldstone bo- 
son(s), while that of the amplitude to the Higgs field, fol- 
lowing the terminology of the electroweak theory. Stim- 
ulated by the role of the phase fluctuation in the un- 
conventional superconducting phase [l47l | in condensed 
matter, we formulat e the 2SC phase in the nonlinear re- 
alization framework jl48j in order to naturally take into 
account the contribution from the phase fluctuation or 
pseudo Nambu-Goldstone current(s). 

In the 2SC phase, the color symmetry G — SU{3)c 
breaks to H ^ SU{2)c- The generators of the residual 
SU{2)c symmetry H arc {S"" = T"} with a 1, 2, 3 and 
the broken generators {X^ = T''+'^} with b = 1, • • • ,5. 
(More precisely, the last broken generator is a combina- 
tion of Tg and the generator 1 of the global U{1) symme- 
try of baryon number conservation, B = [1 + •\/3T8)/3 
of generators of the global U{1)b and local S'[/(3)c sym- 
metry. ) 

The cosct space G/H is parameterized by the group 
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elements 



V{x) = exp 



(53) 



which is the NcNf x iV^jVf- diin ensional Maurer-Cartan 
one- form introduced in Ref. |l48j . The hnear order of the 
Nambu-Goldstone currents and has the exphcit 
form of 



here operator V is unitary, and V^^ = and Lpa{a = 
4, • • • ,7) and ipg, are five Nambu-Goldstone diquarks, and 
we have not consider the topologically nontrivial case and 
therefore V{x) can be expanded uniformly according to 
the powers of tp's. In fact, V{x) is alway topologically 
trivial for a configuration of ip^s that has a finite energy 
because of the trivial homotopy group ■K2{SU{3)/SU (2)). 

Introducing a new quark field x, which is connected 
with the original quark field q in Eq. ()46p through a 
nonlinear transformation, 



and the charge-conjugate fields transform as 

qc = V* xc , qc = xc . 



(54) 



(55) 



The advantage of transforming the quark fields is that 
this preserves the simple structure of the terms coupling 
the quark fields to the diquark sources, 

gcA+9 = xc*^X , q^~qc^X^~Xc- (56) 
In the Nambu-Gor'kov space of the new spinors 



X 



- I X 
XC 



, X=ix,Xc), 



(57) 



the nonlinear realization of the original Lagrangian den- 
sity Eq. takes the form of 
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Here the explicit form of the free propagator for the new 
quark field is 
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and 



if^-/l7o+7M^<^- (61) 



Comparing with the free propagator in the original La- 
grangian density, the free propagator in the non-linear re- 
alization framework naturally takes into account the con- 
tribution from the Nambu-Goldstone currents or phase 
fluctuations, i.e., 
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The advantage of the non-linear realization framework 
Eq. (j58p is that it can naturally take into account 
the contribution from the phase fluctuations or Nambu- 
Goldstone currents. The task left is to find the correct 
ground state by exploring the stability against the fluc- 
tuations of the magnitude and the phases of the order 
parameters. The free energy ri(Vp, /.t, /ig, /ig) can be 
evaluated directly and it takes the form of 
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In the 2SC phase, the color symmetry SU{Z)c is spon- 
taneously broken to SU{2)c and diquark field obtains a 
nonzero expectation value. Without loss of generality, 
one can always assume that diquark condenses in the 
anti-blue direction, i.e., only red and green quarks par- 
ticipate the Cooper pairing, while blue quarks remains 
as free particles. The ground state of the 2SC phase is 
characterized by (A^) = A, and (A^) = 0, (A^) = 0. 

Considering the fluctuation of the order parameter, the 
diquark condensate can be parameterized as 
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(66) 



where $''(a;) = (0, 0, A -I- H{x)) is the diquark field in 
the nonlinear realization framework, i^aT^i are Nambu- 
Goldstone bosons, and H is the Higgs field. 

Expanding the diquark field around the ground 
state: = (0, 0, A), the free-energy of the system takes 
the following expression as 
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There are three contributions to the free-energy, the 
mean-field approximation free-energy part has the 
form of 
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the free-energy from the Higgs field VIh has the form of 
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and the free-energy from the Nambu-Goldstone currents 
flNG has the form of 
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with 



^^(K) = gA^^{K)Ta~V^{K) , 
^(A-) = -gA^(if)Tj-fFc^(i^) 



(72) 
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where the inverse propagator takes the form of 
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(74) 

The quasi-quark propagator at mean-field approximation 
has the form of 



G+ E- 



E+ G- 



(75) 



and its explicit expression of the Namb u- Go rkov compo- 
nents of 5a/ has been derived in Ref. |l20l | . 

The Matsubara self-energy functions IIh{K) of the 
Higgs field and that of the Goldstone fields (obtained af- 
ter the sum over pq and integral over p in Eq. (|7ip ) can be 
continuated to real frequency following the standard pro- 
cedure. Because of the gapless excitations, the values of 
these functions at zero frequency and zero mometum de- 
pends the order of the limit. In this work, we shall restrict 
our attention to static fluctuations only which amounts 
to replace H{K), A{K) and ip{K) of Eqs. §^ and ^ 

by \/TH{k)Skofi, VTA{k)5ko,o and VT(p{k)5koft. Thus 
the long wavelength limit of the self-energy functions dis- 
cussed below corresponds to the limit limj_j_p limfcQ_i.o. 



Sarma instability 

In the mean-field approximation, the fre e-en ergy for 
u, d quarks in /3-equilibrium takes the form jll8j |: 
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which is the same as Eq. (|48| by neglecting the contribu- 
tion of free electrons and taking zero quark mass in the 
color superconducting phase. 

At zero temperature, the mean-field free-energy has 
the expression of 



M 



AGd 27r2 127r2 Utt"^ 
'■^ p^dp 







^ {Vip + f^y- + A2 + v(p - fir + 

(77) 

As we already knew that, with the increase of mis- 
match, the ground state will be in the gapless 2SC phase 
when A < dfi, the thermodynamical potential of which 
is given by 



Ao 



(78) 

where fl^^j is the normal phase thermodynamic potential. 
Aq the solution to the gap equation in the absence of 
mismatch, (5/i = 0. The solution to the gap equation 
reads 



(79) 



Th e gap less phase is in principle a metastable Sarma 
i.e., the free-energy is a local maximum with 



state 
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respect to the gap parameter A. We have 
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(80) 



The weak coupling approximation is employed in deriving 
Eqs. (mi) and ^ from Eq. dTT]), which assumes that Aq, 
A and Sfi are much smaller than fi and A — fi. The same 
approximation will be applied throughoutthe paper. 

Nambu-Goldstone currents generation and the 
LOFF state 

The quadratic action of the Goldstone modes in the 
long wavelength limit can be written down with the aid 
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of the Meissner masses evaluated in Ref. jl20l |. We find 
that 



NG 
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+ higher orders . 
where toi = m2 = rn^ = 0, 
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It was found that at zero temperature, with the increase 
of mismatch, for five gluons witli a = 4,5,6,7,8 corre- 
sponding to broken generator of SU{3)c, their Meissner 
screening mass squares become negative |l20l |. This in- 
dicates the development of the condensation of 
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It can be interpreted as the S|)ontaneous generation of 
Nambu-Goldstone currents < Vip'^ >^ |l39| . or 



gluon condensation X]a=4 < Il40l. I t can also 

be interpreted colored-LOFF state 11411 with the 
plane- wave order parameter 



A{x) = Ae'^'=4 



(85) 



Higgs instability 



We discussed the two known instabilities induced by 
mismatch, i.e., the Sarma instability and chromomag- 
netic instability, respectively. It is found that there is 
another instability, which is related to the Higgs field, 
and we call this instability " Higgs instability" . 

The free-energy from the Higgs field can be evaluated 
and takes the form of 
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Evaluating the one- loop quark-quark bubble H/f(fc), 
we obtain that 
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where the functions /(/cjJ/i) and J{k\6fi) arc given by 
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with Re^ (cj„ + i(5/i)2 + A2 > 0. Here, the summation 
over n is the frequency summation at finite temperature 
field theory, with uj„ = (2n + l)nT. At T=0, the free 
energy of the Higgs field has the form of 
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It follows from Eqs. ^ and 1^ that the Higgs field 
becomes unstable in the gapless phase when dfi > A, Ah 
in the gapless phase is shown in Fig. [19] by the red dash- 
dotte d lin e. The Higgs instability was also considered in 
Ref. |l43| , where it was called as " amplitude instability" . 
It has to be pointed out that we got different expressions 
for the coefficients of the gradient term. For k >> A, we 
have 
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Therefore the Higgs instability disappears for sufficiently 
large momentum. The form factor H^f(fc) for arbitrary 
momentum is plotted using red dashed line in Fig. 1201 for 
a typical value of the mismatch parameter A/Sfi = 1/2. 
We notice that the Higgs instability becomes stronger for 
nonzero momentum. 

Charge neutrality condition 

Unless a competing mechanism that results in a pos- 
itive contribution to the inhomogeneous blocks of the 
stability matrix. The Higgs instability will prevent gap- 
less superfluidity/ superconductivity from being imple- 
mented in nature. In the system of imbalanced neutral 
atoms, such a mechanism is not likely to exist and this 
contributes to the reason why the BP state has never 
been observed there. For the quark matter being consid- 
ered, however, the positive Coulomb energy induced by 
the Higgs field of electrically charged diquark pairs has 
to be examined. 

The second order derivative of the Hemlholtz free en- 
ergy for g2SC reads: 
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FIG. 19: Squared values of the gluon Debye (upper panel) 
and Meissner (lower panel) screening masses, devided by 
rrig = Aasp.^ /Sn, as functions of the dimensionless parameter 
A/Sfx. The red solid line denotes the results for the giuons 
with A = 1,2, 3, the green long-dashed line denotes the results 
for the giuons with A = 4,5, 6, 7, and the blue short-dashed 
line denotes the r esult s for the gluon with A — 8. This figure 
is taken from Ref. Il42ll . 



and the stability implemented by the charge neutrality 
implies that 



f— ) 



> 0. 



(95) 



While the global charge neutrality is maintained, an 
inhomogeneous (j) will induce local charge distribution. 
The corresponding Coulomb energy 



1 ^ 6p{k)*5p{k) 

fc#0 



(96) 



should be considered where Sp(k) is the Fourier com- 
ponent of the (/)-induced charge density and to^ is the 
Coulomb polarization function ( Debye mass at fc = ). 
We have 

ml{k)^-'-^Yl ^^^oQSiP + K)^oQSiP) (97) 
^ p 



and 



Sp{k) = K{k)H{k), 



(98) 



where 



^(^) = ^ E trToQ5(P + K)j5e^p2SiK) (99) 



with K = (0, k) and Q the electric charge operator. We 
have 



Q = P3ia + bTs) 



(100) 



with a = 1/6 and 6=1/2 for the quark matter consisting 
of u and d flavors. 

The modified Higgs self-energy has the form of 



flnik) 



dH*[k)dH(ky ri^ 

d^n \ K*{k)K{k) 



- + 



dH*{k)dH{k)^ f^. ' k^ + mlik) 

K*{k)K{k) 

= ^"^^^ + k^ + miiky 



(101) 



The stability of the system with respect to the Higgs field 
requires that nff(fc) > for all k. 

The form factor K{q) and the momentum dependent 
Debye mass square can be calculated explicitly at weak 
coupling, i.e. Sfi << p and k << p. We find that 



(102) 



with 



Kik\6n) = iATj2 



dx 



1 (w„ + i(5//)2 + A2 ifc 



TTi7i(103) 



and 



m^{k) 



{k\6p) (104) 



with the first term the Debye mass of the normal phase 
and I{q\Sp) the function defined in the section III. 

Sarma instability at zero momentum limit can 
be removed by Coulomb energy: 



Notice that 

Imi lim — 

k^ov^oo \dH*(k)dH{k)KQ 



and 



We have 



lim 



,(0) = 



k(0) = e 



/drie\ 



k^o ^ dH* {k)dH{k) ^ t^e V dlS? J ' 



( 



(106) 



(107) 



(108) 



and the charge neutrality stabilize also the inhomoge- 
neous Higgs field with the momentum much smaller 



20 



than the inverse coherence length and the inverse Debye 
length. 

In the static long- wave length limit, we have 



I'D 



(0) = 



1 



and 



(109) 



(110) 



It follows from Eqs. ([801), (fTOT|) . (fT09)) and (fTTO)) that the 
Higgs self-energy including Coulomb energy correction 



iff =IIh(0) 



4(62 _ 3a2)^2(^^ _ ^5^2 _ a2) 



7r2 [3a2 ^<5Ai2 - A2 + fe2 (2^^i + ^5fi^ - A^)] 

(111) 

is always positive for the whole range of g2SC state. It 
means that the Sarma instability in the gapless phase can 
be cured by Coulomb energy. This is shown in Fig. I19[ 
where the red dash-dotted line indicates the Ajj and the 
green solid line indicates Ah- 
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FIG. 20: The function Unik) (red dashed line), n_tf (fc) (black 
solid line) and the Coulomb energy (green dash-dotted line) 
as functions of scaled-momentum k/A, in the case of Sp, = 2 A 
and (e^^t^)/(47rA^) = 1. This figure is taken from Ref [li^l. 



Higgs instability at nonzero momentum cannot 
be removed by Coulomb energy: 

While the Sarma instability in g2SC phase can be 
cured by Coulomb energy under the constraint of charge 
neutrality condition, it is not sufficient for the system 
to be stable, even if the chromomagnetic instabilities are 
removed, say by gluon condensation. One has to explore 
the Higgs instability by calculating the self-energy func- 
tion n(fc) in the whole momentum space, which amounts 
to value the three basic functions /(fc|(5/i), J(fc|(5/i) and 
K{k\5^i) defined in ^ and (ITU51) . 

Fig. [20] shows the Higgs self-energy linik) (the red 
dashed line) the Coulomb corrected Higgs self-energy 
n//(fc) (the black solid line) and the Coulomb en- 
ergy Ecoui (the green dash-dotted line) as functions of 
scaled- momentum k/A, in the case of J/i = 2 A and 
(e2/l2)/(47rA2) 1. 

Eventhough the Higgs instability can be removed by 
the Coulomb energy for small momenta, it returns for 
intermediate momenta. This phenomenon persists for a 
wide range of gap magnitude, < A < 0.866(5/i and 
for all strength of the Coulomb interaction, measured by 

-2 

the dimensionless ration ?/ = . Within the narrow 
range 0.8665/.t < A < Sfi, the Higgs instability could 
be removed if the Coulomb interaction were sufficiently 
strong. For 77 = 1, We found that linik) > for all k if 
0.998(5^ < A < Sfi. In terms of the values of the param- 
eters of NJL model, we have aefl^ < Sjj,. Therefore the 
electric Coulomb energy cannot cure the Higgs instability 
for a realistic two flavor quark matter. 



Negative Hff (fc) indicates the Higgs mode is unstable 
and will decay jl49l |. It is noticed that Unik) reaches its 
minimum at a momentum, i.e., k ~ 4A, which indicates 
that a stable state may develop around this minimum, we 
characterize this momentum as fcmm • The inverse k~]„ 
is the typical wavelength for the unstable mode jl50l |. 
If mixed phase can be formed, the typical size I of the 
2SC bubbles should be as great as k^]^, i.e., / ~ fc~ 
[l50j | , which turns out to becomparable to the coherence 
length of 2SC in accordance with Eq. ([75)) . (In the case 
of homogeneous superconducting phase, kmin = 0, and 
Z — )- 00.) Considering that the coherence length ^ of a 
superconductor is proportianl to the inverse of the gap 
magnitude, i.e., ^ ~ A~^, therefore, a rather large ratio 
of kjnin/ A means a rather small ratio of l/£,. When l/^ < 
1. a phase separation state is more favorable. 

The reader has to keep in mind that our result of the 
Higgs instability only indicates some kind of inhomoge- 
ncous states whose typical length scale is comparable to 
the coherence length of 2SC. Further insight on the struc- 
ture of the inhomogeneity cannot be gained without ex- 
ploring the higher order terms of the nonlinear realiza- 
tion ([66]) . A more direct approach to obtain the favorite 
structure of the ground state is to compare the free en- 
ergy of various candidate states, which include the mixed 
phase, the single-plane wave FF state, striped LO state 
and multi-plane wave states. The Coulomb energy and 
the gradient energy have to be estimated reliably. We 
leave this analysis as a future project. 

In the system of imbalanced neutral atoms, the Higgs 
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instability persists and induces spatial non-uniform phase 
separation state. This explains why imbalanced cold 
atom experiments did not observe LOFF state rather 
showed strong evidence of phase separation. For the 
2-flavor quark matter being considered, the electric 
Coulomb interaction is not strong enough to compete 
with the Higgs instability. 

IV. CONCLUSION 

I have introduced several topics of QCD phase struc- 
ture at high temperature and high density: the properties 
of strongly interacting quark gluon plasma, searching for 
the critical end point and the gapless color superconduc- 
tor. 

I give a brief introduction on the discovery of sQGP. It 
has been believed for more than 30 years that the QGP 
created at heavy-ion collisions should be weakly inter- 
acting gas system. However, a very small shear viscosity 
over entropy density ratio 77/s is required to fit the RHIC 
data of elliptic flow V2 ■ It is in contrary to the large value 
of rj/s given by perturbative QCD calculation. This is 
a strong evidence that the deconfined matter created at 
RHIC is strongly interacted. The AdS/CFT duality gives 
a lower bound rj/s = I/Att. Therefore, it is conjectured 
that the sQGP created at RHIC might be the most per- 
fect fluid observed in nature. However, a perfect fluid 
should have both vanishing shear and bulk viscosities. 
The property of the bulk viscosity of sQGP need to be 
studied before one draws the flnal conclusion. 

Recent studies show that the bulk viscosity over en- 
tropy density ratio ^/s rises up near phase transitions. 
The result from QCD effective models show that C/s be- 



haves differently for different orders of phase transitions: 
for Ist-order phase transition, it rises sharply and show 
a divergent behavior, for the 2nd-order phase transition, 
it shows an upward cusp at Tc, for the case of crossover, 
the cusp becomes smooth. It is discussed that the sharp 
rising bulk viscosity will lead to the breakdown of hydro- 
dynamics and affect the hadronization. Therefore the 
critical end point might be located through the observ- 
ables which arc sensitive to the ratio of bulk viscosity 
over entropy density. 

The status of gapless color superconductor is reported. 
The chromomagnetic instability, the Sarma instability 
and Higgs instability are clarified. In the gapless color su- 
perconducting phase, both the phase part and magnitude 
part of the order parameter will develop instabilities: The 
phase part develops into the chromomagnetic instability, 
which induces the plane- wave state; The magnitude part 
develops the Sarma instability and Higgs instability, the 
Sarma instability can be competed with charge neutral- 
ity condition, while the Higgs instability cannot be cured 
by Coulomb interaction, and induces the inhomogeneous 
state. 

V. ACKNOWLEDGEMENTS 



I thank the collaboration with J.W.Chen, B.C.Li, 
H.Mao, D.L.Yang for the sQGP part, and I.Giannakis, 
D.Hou, H.C.Ren and I. Shovkovy for the gapless color 
superconductor part. The work is supported by CAS 
program " Outstanding young scientists abroad brought- 
in", CAS key project KJCX3-SYW-N2, NSFC10735040, 
NSFC10875134, and the K.C.Wong Education Founda- 
tion, Hong Kong. 



[1] D.J. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973) 

1343; H.D. Politzer, Phys. R ev. Lett. 30 (1973) 1346. 
[2] E. Klempt, |hep-ph/0404270 ]j 

[3] T.D. Lee and G.C. Wick, Phys. Rev. D 9, 2291 (1974). 
[4] J.C. Collins and M.J. Perry, Phys. Rev. Lett. 34, 1353 
(1975). 

[5] G. Baym and S.A. Chin, Phys. Lett. B 62, 241 (1976). 

[6] E.V. Shuryak, Phys. Lett. B 78, 150 (1978). 

[7] F. Karsch and E. Laermann, Phys. Rev. D 50, 6954 

(1994); F. Karsch, Nucl. Phys. A 698, 199 (2002). 
[8] E.V. Shuryak and I. Zalied, Phys. Rev. C 70, 021901 

(2004); E. Shuryak, Prog. Part. Nucl. Phys. 53, 273 

(2004) ; E.V. Shuryak and I. Zahed, Phys. Rev. D 70, 
054507 (2004); E.V. Shuryak, Nucl. Phys. A 750, 64 

(2005) . 

[9] D.H. Rischke, Prog. Part. Nucl. Phys. 52, 197 (2004). 
[10] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 
(2005). 

[11] P. Jacobs and X.N. Wang, Prog. Part. Nucl. Phys. 54, 

443 (2005). 
[12] U.W. Heinz, [hep^h /0407360' 

[13] H.R. Jaqaman, A.Z. Mekjian, and L. Zamick, Phys. 



Rev. C 29, 2067(1984); R.K. Su, F.M. Lin Phys. Rev. 
C 39, 2438(1989). 

[14] U.W. Heinz, P. R. Subramanian, H. Stocker and 
W. Greiner, J. Phys. G 12, 1237 (1986). 

[15] S.C. Frautschi, Asymptotic freedom and color super- 
conductivity in dense quark matter, Proc.of the Work- 
shop on Hadronic Matter at Extreme Energy Density, 
N. Cabibbo (Editor), Erice, Italy (1978). 

[16] F. Barrois, Nucl. Phys. B129 (1977), 390. 

[17] J. Bardeen, L.N. Cooper, J.R. Schrieffer Phys. Rev. 
106:162 (1957); Phys. Rev. 108:1175 (1957). 

[18] D. Baihn and A. Love, Phys. Rep. 107, 325(1984). 

[19] R. Rapp, T. Schafer,E.V. Shuryak and M. Velkovsky, 
Phys.Rev.Lett. 81,53(1998); M. Alford, K. Rajagopal 
and F. Wilczek, Phys.Lett.B 422,247( 1998). 

[20] K. Rajagopal and F. Wilczek, hep-ph/0011333j 
D.K. Hong, Acta Phys. Polon. B 32, 1253 (2001); M. Al- 
ford, Ann. Rev. Nucl. Part. Sci. 51, 131 (2001); G. Nar- 
duUi, Riv. Nuov o Cim. 25N3, 1 (2002); T. Schafer, 
hep-ph/0304281[ M. Bu balla, Phys. Rept. 407, 205 
(2005); H.C. Ren, 'hep-ph/0404074 M. Huang, Int. J. 
Mod. Phys. E 14, 675 (2005); I. A. Shovkovy, Found. 



22 



Phys. 35, 1309 (2005); M. G. Alford, A. Schmitt, K. Ra- 
jagopal and T. Schafer, Rev. M od. Phys. 80, 1455 

(2008) ; Q. Wang. EXi^:0912.2485l [nucl-th]. 

[21] P. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0305, 
051 (2003). 

[22] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. 
Lett. 86^ 4783 (2001), P. Huovinen, P. F. Kolb, 
U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, 
Phys. Lett. B 503, 58 (2001), T. Hhano, U. W. Hemz, 
D. Kharzeev, R. Lacey and Y. Nara, Phys. Lett. B 
636, 299 (2006), P. Romatschke and U. Romatschke, 
Phys. Rev. Lett. 99, 172301 (2007), H. Song and 
U. W. Heinz, Phys. Lett. B 658, 279 (2008), H. Song 
and U. W. Heinz, Phys. Rev. C 77, 064901 (2008). 
D. Teaney, Phys. Rev. C 68, 034913 (2003). 
A. Nakamura and S. Sakai, Phys. Rev. Lett. 94, 072305 
(2005). 

L Arsene et al. [BRAHMS Collaboration], Nucl. Phys. 
A 757, 1 (2005), K. Adcox et al. [PHENIX Collabora- 
tion], Nucl. Phys. A 757, 184 (2005), B. B. Back et a/., 
Nucl. Phys. A 757, 28 (2005), J. Adams et al. [STAR 
Collaboration], Nucl. Phys. A 757, 102 (2005). 
M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 
(2005). 

G. Policastro, D. T. Son and A. O. Starinets, Phys. 
Rev. Lett. 87, 081601 (2001), P. Kovtun, D. T. Son and 

A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005). 
P. Arnold, C. Dogan and G. D. Moore, Phys. Rev. D 
74, 085021 (2006). 

D. Kharzeev and K. Tuchin, JHEP 0809, 093 (2008). 

F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B 
663, 217 (2008). 

H. B. Meyer, Phys. Rev. Lett. 100, 162001 (2008). 
K. Huebner, F. Karsch and C. Pica, Phys. Rev. D 78, 
094501 (2008). 

K. Paech and S. Pratt, Phys. Rev. C 74, 014901 (2006). 

B. C. Li and M. Huang, Phys. Rev. D 78, 117503 (2008); 

B. C. Li and M. Huang, Phys. Rev. D 80, 034023 (2009). 
J. W. Chen and J. Wang, Phys. Rev. C 79, 044913 

(2009) . 

C. Sasaki and K. Redlich. larXiv:0811.4708 [hep-ph]. 

G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Leg- 
eland, M. Lutgemeier and B. Petersson, Nucl. Phys. B 
469, 419 (1996). 

M. Cheng et al, Phys. Rev. D 77, 014511 (2008). 

D. Fernandez-Fraile and A. G. Nicola, Phys. Rev. 
Lett. 102, 121601 (2009); D. Fernandez-Fraile and 
A. G. Nicola, Eur. Phys. J. C 62, 37 (2009). 
G. Torrieri, B. Tomasik and L Mishustin, Phys. Rev. C 
77, 034903 (2008); G. Torrieri and L Mishustin, Phys. 
Rev. C 78, 021901 (2008). 

R. J. Fries, B. MuUer and A. Schafer, Phys. Rev. C 78, 
034913 (2008). 

T. Brasoveanu, D. Kharzeev and M. Martinez, Lect. 
Notes Phys. 785, 341 (2010). 
J. I. Kapusta, arXiv:0809.3746 [nucl-th]. 
Z. Fodor and S. D. Katz, J. High Energy Phys. 0203, 
014 (2002), P. de Forcrand and O. Philipsen, Nucl. Phys. 
Proc. Suppl. 129, 521 (2004), F. Karsch, C. R. AUton, 
S. Ejiri, S. J. Hands, O. Kaczmarek, E. Laermann, and 
C. Schmidt, Nucl. Phys. Proc. Suppl. 129, 614 (2004). 
[45] M. Asakawa and K. Yazaki, Nucl. Phys. A 504, 668 
(1989), A. Barducci, R. Casalbuoni, G. Pettini and 
R. Gatto, Phys. Rev. D 49, 426 (1994), J. Berges and 



K. Rajagopal, Nucl. Phys. B 538, 215 (1999), M. A. Ha- 
lasz, A. D. Jackson, R. E. Shrock, M. A. Stephanov and 
J. J. M. Verbaarschot, Phys. Rev. D 58, 096007 (1998), 
P. Zhuang, M. Huang and Z. Yang, Phys. Rev. C 62, 
054901 (2000), O. Scavenius, A. Mocsy, L N. Mishustin 
and D. H. Rischke, Phys. Rev. C 64, 045202 (2001), 
N. G. Antoniou and A. S. Kapoyannis, Phys. Lett. B 
563, 165 (2003). 
[46] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 
(1984). 

[47] Y. Hatta and T. Ikeda, Phys. Rev. D 67, 014028 (2003); 
D. T. Son and M. A. Stephanov, Phys. Rev. D 70, 
056001 (2004); H. Fujii, Phys. Rev. D 67, 094018 (2003); 
H. Fujh and M. Ohtani, Phys. Rev. D 70, 014016 (2004); 
H. Fuju and M. Ohtani, Prog. Theor. Phys. Suppl. 153, 
157 (2004); K. Yagi, T. Hatsuda, Y. Miake, "Quark- 
Gluon Plasma, from big bang to little bang" , Cambridge 
University Press. 

[48] M. A. Stephanov, K. Rajagopal, and E. V. Shuryak, 
Phys. Rev. Lett. 81, 4816 (1998). M. A. Stephanov, 
K. Rajagopal, and E. V. Shuryak, Phys. Rev. D 60, 
114028 (1999). 

[49] R. A. Lacey et al, Phys. Rev. Lett. 98, 092301 (2007). 

[50] L. P. Csernai, J. L Kapusta and L. D. McLerran, Phys. 
Rev. Lett. 97, 152303 (2006). 

[51] S. S. Gubser, A. Nellore, S. S. Pufu and F. D. Rocha, 
Phys. Rev. Lett. 101, 131601 (2008). S. S. Gubser, 
S. S. Pufu and F. D. Rocha, JHEP 0808, 085 (2008). 

[52] J. W. Chen, M. Huang, Y. H. Li, E. Nakano and 
D. L. Yang, Phys. Lett. B 670, 18 (2008)^ 

[53] H. Mao, J. Jin and M. Huang, arXiv:0906 T324 | [hep-ph] . 
to appear in J. Phys. G. 

[54] J.M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. 
D 10, 2428 (1974). 

[55] S. Jeon, Phys. Rev. D 52, 3591 (1995); S. Jeon and L. 
Yaffe, Phys. Rev. D 53, 5799 (1996)^ 

[56] B. J. Schaefer and M. Wagner, larXiv:0812. 28551 [hep- 
ph]. 

[57] K. Fukushima, Phys. Lett. B 591, 277 (2004). 

[58] C. Ratti, M. A. Thaler and W. Weise, Phys. Rev. D 73, 

014019 (2006). 
[59] Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004). 
[60] W. j. Fu, Z. Zhang and Y. x. Liu, Phys. Rev. D 77, 

014006 (2008). 

[61] M. Ciminale, R. Gatto, N. D. Ippolito, G. NarduUi and 
M. Ruggieri, Phys. Rev. D 77, 054023 (2008). 

[62] L. McLerran and R. D. Pisarski, Nucl. Phys. A 796, 83 
(2007). 

[63] G. D. Moore and O. Saremi, JHEP 0809, 015 (2008). 
[64] A. Onuki, Phys. Rev. E 55, 403 (1997). 
[65] L.N. Cooper, Phys. Rev. 104, 1189 (1956). 
[66] W. Meissner and R. Ochsenfeld, Naturwiss. 21, 787 
(1933). 

[67] D.T. Son, Phys. Rev. D 59, 094019 (1999). 
[68] T. Schafer and F. Wilczek, Phys. Rev. D 60, 114033 
(1999). 

[69] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 
074017 (2000). 

[70] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 

051501 (2000). 
[71] D.K. Hong, Phys. Lett. B 473, 118 (2000). 
[72] D.K. Hong, V.A. Miransky, LA. Shovkovy and 

L.C.R. Wijewardhana, Phys. Rev. D 61, 056001 (2000); 

[Erratum- ibid. D 62, 059903 (2000)]. 



23 



[73] W.E. Brown, J.T. Liu and H.C. Ren, Phys. Rev. D 61, 
114012 (2000); W.E. Brown, J.T. Liu and H.C. Ren, 
Phys. Rev. D 62, 054016 (2000); W.E. Brown, J.T. Liu 
and H.C. Ren, Phys. Rev. D 62, 054013 (2000). 
[74] C. Manuel, Phys. Rev. D 62, 114008 (2000). 
[75] R. Rapp, T. Schafer, E.V. Shuryak and M. Velkovsky, 
Annals Phys. 280, 35 (2000). 
6] J. Berges and K. Rajagopal, Nucl. Phys. B 538, 215 
(1999). 

7] G.W. Carter and D. Diakonov, Phys. Rev. D 60, 016004 
(1999). 

S] T.M. Schwarz, S.P. Klevansky and G. Papp, Phys. Rev. 

C 60, 055205 (1999). 
9] R. Nebauer and J. Aichehn, Phys. Rev. C 65, 045204 
(2002) 

3] M. Buballa, J. Hosek and M. Oertel, Phys. Rev. D 65, 
014018 (2002); M. BubaUa and M. Oertel, Nucl. Phys. 
A 703, 770 (2002). 
1] D. Ebert, K.G. Klimenko and H. Toki, Phys. Rev. 
D 64, 014038 (2001); D. Ebert, V.V. Khudyakov, 
V.C. Zhukovsky and K.G. Klimenko, Phys. Rev. D 65, 
054024 (2002). 

2] M. Huang, P.P. Zhuang and W.Q. Chao, Phys. Rev. D 

65, 076012 (2002). 
3] T. Schafer and F. Wilczek, Phys. Lett. B 450, 325 

(1999). 

4] N.J. Evans, S.D.H. Hsu and M. Schwetz, Nucl. Phys. 
B 551, 275 (1999); N.J. Evans, S.D.H. Hsu and 
M. Schwetz, Phys. Lett. B 449, 281 (1999); S.D.H. Hsu 
and M. Schwetz, Nucl. Phys. B 572, 211 (2000). 
[85] D.H. Rischke, Phys. Rev. D 62, 034007 (2000); 
D.H. Rischke and LA. Shovkovy, Phys. Rev. D 66, 
054019 (2002). 

[86] M.G. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. 

B537, 443 (1999). 
[87] T. Schafer, Nucl. Phys. B 575, 269 (2000). 
[88] N.J. Evans, J. Hormuzdiar, S.D.H. Hsu and M. Schwetz, 

Nucl. Phys. B 581, 391 (2000). 
[89] LA. Shovkovy and L.C.R. Wijewardhana, Phys. Lett. B 

470, 189 (1999). 
[90] D.H. Rischke, Phys. Rev. D 62, 054017 (2000). 
[91] R. Casalbuoni and R. Gatto, Phys. Lett. B 464, 111 

(1999) . 

[92] D. T. Son and M. A. Stephanov, Phys. Rev. D 61, 
074012 (2000); [Erratum-ibid. Phys. Rev. D 62, 059902 

(2000) ]. 

[93] T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 
(1999). 

[94] T. Schafer, Phys. Rev. D 62, 094007 (2000). 

[95] M. Buballa, J. Hosek and M. Oertel, Phys. Rev. Lett. 

90, 182002 (2003). 
[96] A. Schmitt, Q. Wang and D.H. Rischke, Phys. Rev. D 

66, 114010 (2002). 

[97] A. Schmitt, Q. Wang and D.H. Rischke, Phys. Rev. 
Lett. 91, 242301 (2003); A. Schmitt, Q. Wang and 
D.H. Rischke, Phys. Rev. D 69, 094017 (2004). 

[98] A. Schmitt, nucl-th/0405076 

[99] A.L Larkin, YuN. Ovchinnikov, Zh. Eksp. Tear. Fiz. 47 
(1964), 1136; translation: Sov. Phys. JETP 20 (1962), 
762. 

[100] P. Fulde, R.A. FerreU, Phys. Rev. 135, A550 (1964). 
[101] M.G. Alford, J.A. Bowers and K. Rajagopal, Phys. Rev. 

D 63, 074016 (2001). 
[102] J.A. Bowers, J. Kundu, K. Rajagopal and E. Shuster, 



Phys. Rev. D 64, 014024 (2001). 
[103] A.K. Leibovich, K. Rajagopal and E. Shuster, Phys. 

Rev. D 64, 094005 (2001). 
[104] J. Kundu and K. Rajagopal, Phys. Rev. D 65, 094022 

(2002). 

[105] J.A. Bowers and K. Rajagopal, Phys. Rev. D 66, 065002 
(2002). 

[106] J.A. Bowers, hep-ph/0305301 

[107] I. Giannakis, J.T. Liu and H.C. Ren, Phys. Rev. D 66, 
031501 (2002). 

[108] R. Casalbuoni and G. Nardulh, Rev. Mod. Phys. 76, 
263 (2004); R. Casalbuoni, M. Ciminale, M. MannareUi, 
G. NarduUi, M. Ruggieri and R. Gatto, Phys. Rev. 
D 70, 054004 (2004); R. Casalbuoni, R. Gatto, 
M. MannareUi, G. NarduUi and M. Ruggieri, Phys. Lett. 
B 600, 48 (2004). 

[109] T. Schafer, Phys. Rev. Lett. 85, 5531 (2000); P.P. Be- 
daque and T. Schafer, Nucl. Phys. A 697, 802 (2002). 

[110] D.B. Kaplan and S. Reddy, Phys. Rev. D 65, 054042 

(2002) . 

[Ill] T. Schafer, D.T. Son, M.A. Stephanov, D. Toublan and 
J.J.M. Verbaarschot, Phys. Lett. B 522, 67 (2001). 

[112] V.A. Miransky and LA. Shovkovy, Phys. Rev. Lett. 88, 
111601 (2002). 

[113] M. Rho, A. Wirzba and L Zahed, Phys. Lett. B 473, 
126 (2000); 

[114] D.K. Hong, T. Lee and D.P. Min, Phys. Lett. B 477, 
137 (2000). 

[115] C. Manuel and M.H.G. Tytgat, Phys. Lett. B 479, 190 

(2000) . 

[116] T. Schafer, Phys. Rev. D 65, 074006 (2002). 

[117] A. Kryjevski, D.B. Kaplan and T. Schafer, Phys. Rev. 

D 71, 034004 (2005). 
[118] I. Shovkovy and M. Huang, Phys. Lett. B 564, 205 

(2003) . 

[119] M. Huang and I. Shovkovy, Nucl. Phys. A 729, 835 

(2003) . 

[120] M. Huang and LA. Shovkovy, Phys. Rev. D 70, 051501 

(2004) ; M. Huang and LA. Shovkovy, Phys. Rev. D 70, 
094030 (2004). 

[121] M. Alford, C. Kouvaris and K. Rajagopal, Phys. Rev. 

Lett. 92, 222001 (2004); M. Alford, C. Kouvaris and 

K. Rajagopal, Phys. Rev. D 71, 054009 (2005). 
[122] K. lida, T. Matsuura, M. Tachibana and T. Hatsuda, 

Phys. Rev. Lett. 93, 132001 (2004). 
[123] S.B. Ruster, LA. Shovkovy and D.H. Rischke, Nucl. 

Phys. A 743, 127 (2004). 
[124] K. Fukushima, C. Kouvaris and K. Rajagopal, Phys. 

Rev. D 71, 034002 (2005). 
[125] A. Kryjevski and T. Schaefer, Phys. Lett. B 606, 52 

(2005) . 

[126] M. Alford and K. Rajagopal, JHEP 0206, 031 (2002). 
[127] A.W. Steiner, S. Reddy and M. Prakash, Phys. Rev. D 

66, 094007 (2002). 

[128] M. Huang, P.P. Zhuang and W.Q. Chao, Phys. Rev. D 

67, 065015 (2003). 

[129] S.B. Ruster, Diploma thesis, J. W. Goethe-University, 
2003. S.B. Ruster and D.H. Rischke, Phys. Rev. D 69, 
045011 (2004). 

[130] A. Mishra and H. Mishra, Phys. Rev. D 69, 014014 
(2004). 

[131] K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86, 3492 

(2001) . 

[132] G. Sarma, J. Phys. Chem. Solids 24, 1029(1963). 



24 



[133] M.G. Alford, J. Berges and K. Rajagopal, Phys. Rev. 

Lett. 84, 598 (2000). 
[134] W.V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 

(2003); W.V. Liu, F. Wilczek and P. Zoller, Phys. Rev. 

A 70, 033603 (2004); M.M. Forbes, E. Gubankova, 

W.V. Liu and F. Wilczek, Phys. Rev. Lett. 94, 017001 

(2005) . 

[135] E. Gubankova, W.V. Liu and F. Wilczek, Phys. Rev. 
Lett. 91, 032001 (2003). 

[136] E. Gubankova, E. G. Mishchenko and F. Wilczek, Phys. 
Rev. Lett. 94, 110402 (2005). 

[137] S.B. Nam, Phys. Rev. 156, 470 (1967). 

[138] A.K. Geim, S.V. Dubonos, J.G.S. Lok, M. Henini, 
and J.C. Maan, Nature 396, 144 (1998); M. Sigrist, 
T.M. Rice, Rev. Mod. Phys. 67, 503 (1995); Mai Suan 
Li, Phys. Rept. 376, 133(2003). 

[139] D. K. Hong,jhep-ph/0506097; Mci Huang, Int. J. Mod. 
Phys. A 21, 9TEr(2006); Mei Huang, Phys. Rev. D 73, 
045007 (2006); A. Kryjevski, Phys. Rev. D 77, 014018 
(2008); T. Schafer, Phys. Rev. Lett. 96, 012305 (2006); 
A. Gerhold and T. Schafer, Phys. Rev. D 73, 125022 

(2006) . 

[140] E. V. Gorbar, M. Hashimoto and V. A. Miransky, Phys. 
Lett. B 632, 305 (2006); M. Hashimoto, Phys. Lett. B 
642, 93 (2006). 

[141] I. Giannakis and H. C. Ren, Phys. Lett. B 611, 137 
(2005); I. Giannakis and H. C. Ren, Nucl. Phys. B 723, 
255 (2005); I. Giannakis, D. f. Hon and H. C. Ren, Phys. 
Lett. B 631, 16 (2005); 



[142] L Giannakis, D. Hon, M. Huang and H. c. Ren, Phys. 
Rev. D 75, 011501 (2007); L Giannakis, D. Hon, 
M. Huang and H. c. Ren, Phys. Rev. D 75, 014015 
(2007). 

[143] K. lida and K. Fukushima, Phys. Rev. D 74, 074020 
(2006). 

[144] H. Muther and A. Sedrakian, Phys. Rev. D 67, 085024 
(2003) 

[145] E. Nakano, T. Maruyama and T. Tatsumi, Phys. Rev. 

D 68, 105001^ (2003); T. Tatsumi, T. Maruyama and 

E. Nakano, hep-ph/0312351 
[146] A. A. Abrikosov, Sov.Phys.JETP 5, 1174(1957) [Zh. 

Eksp. Theo. Fiz. 32, 1442 (1957)]. 
[147] H. Kleinert, cond-mat/95 03030 M.Franz, Z. Tesanovic, 

Phys.Rev.Lett.87, 257003(2001); M.Franz, Z. 

Tesanovic, O.Vafek, Ph ys.Rev.B 66(2002), 054535; 

A. Mehkyan, Z. Tesanovic, |cond-m at /0408344 
[148] R. Casalbuoni, Z. y. Duan and F. Sannino, Phys. Rev. 

D 62, 094004 (2000); V. A. Miransky, L A. Shovkovy 

and L. C. R. Wijewardhana, Phys. Rev. D 64, 096002 

(2001); D. H. Rischke and L A. Shovkovy, Phys. Rev. 

D 66, 054019 (2002). 
[149] S. R. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 

(1973); S. R. Coleman, Phys. Rev. D 15, 2929 (1977), 

[Erratum-ibid. D 16, 1248 (1977)]. 
[150] E. J. Weinberg and A. q. Wu, Phys. Rev. D 36, 2474 

(1987). 



